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CHAPTER  I 


INTRODUCTION 


History  and  Statement  of  the  Problem 

In  [3],  a  method  was  derived  for  the  estimation  of  the  zeros  of 
a  polynomial  with  real  coefficients,  g(x)  =  ooc    +  f  (x) ,  from  the 
known  zeros  of  f(x),  assuming  that  ck  is  a  small  real  number.   In  [11] 
J.  E.  Maxfield  et  al . ,  established  an  iterative  process  which  con- 
verges to  the  zeros  of  g(x)  and  can  be  used  to  give  a  better  estimate 
of  these  zeros  than  the  previous  method.  The  present  paper  now  aims 
to  present  an  extension  of  [11]: 

Let  (1)  t(x)  be  a  polynomial  with  known  zeros  r  ,  ...  r^. 


n+1 


i-1 


From  t(x)    =     ^     a.x         with  a     =  1   and  the  zeros,    r  ,      of  t(x)    we  wish 
11  n  1 

to  establish   an  iterative  procedure  which  converges  to  the   zeros, 


^r    ^2'    •••    Vp    °^ 

(2)  g(x)    =  a^^^  x""*"     +  f(x)     vrfiere 

n+1  .    , 

(3)  f(x)    =     E  (a.+cy.)x^-^  a^^^=  1,    and 

cy,  .cy^,....   Q'     ^  are  of   small 
1     2'  n+2  .,     , 

magnitude . 

In  this  paper,  the  coefficients  of  f(x),  g(x)  and  t(x),  as 
well  as  the  increments  »  's,  need  not  be  restricted  to  real  numbers. 

Moreover,  as  a  process  for  determining  the  zeros,  r^,  ...  r^, 
of  t(x),  Lin's  Penultimate  Remaindering  Iteration  is  studied,  and  it 


shall  be  shown  that  a  translation  of  the  variable  x  can  control  its 
convergence. 

Application  of  the  Paper 

In  a  problem  for  determining  the  zeros  of  f(x)  with  real  coef- 
ficients, most  of  the  established  methods  that  can  be  applied,  such  as 
Lin's  Penultimate  Remaindering,  are  contaminated  factorization  schemes 
in  the  sense  that  the  degree  of  accuracy  of  one  factor  affects  that  of 
the  succeeding  factors.  For  instance,  in  Penultimate  Remaindering,  the 
depression  of  the  polynomial  after  a  factor  is  obtained  causes  error 
due  to  the  rounding  in  the  division  process  itself,  as  well  as  due  to 
the  inaccuracy,  perhaps,  of  the  first  factor  obtained.  Thus  most  of 
these  methods  give  only  the  estimate  r.  to  the  actual  zeros  r. .  If 
t(x)  is  the  polynomial  having  the  r.s  for  its  true  zeros,  then  each 
coefficient  of  t(x)  differs  from  the  corresponding  coefficient  of  f(x) 
by  some  small  number  a         Using  these  c. 's  and  r.  's  as  input,  the  iter- 
ative procedure  we  derive  here  may  be  applied  as  a  corrective  measure 
to  obtain  the  actual  zeros  r. 's,  of  g(x) |    ^q   which  is  precisely 
f(x). 

Aside  from  being  merely  corrective,  the  said  iteration  can 
also  be  used  as  a  method  by  itself  of  determining  the  zeros  of  a  poly- 
nomial, independent  of  any  other  method.  In  this  case  we  take  for  the 
r.'s  rough  estimates  of  the  zeros  being  sought.  For,  even  though  the 
iteration  shall  be  proved  to  converge  only  for  the  increments  a   ' s 
small  enough  in  absolute  value,  like  other  iterations,  experience  shows 
that  surprisingly,  even  with  very  poor  estimates  of  the  zeros  for  a 
start,  the  iteration  frequently  converges. 


Furthermore,  suppose  the  zeros  of  g(x)  =  a        x         +  f(x), 
where  o"  ^  is  quite  small  in  magnitude,  are  to  be  determined.  We  may- 
truncate  it  temporarily  by  neglecting  the  leading  term  to  obtain  the 
polynomial  f(x)  of  lower  degree.  If  we  can  factor  this  poljTiomial  and 
with  these  zeros  as  r, ' s,  we  may  then  employ  our  iterative  scheme  to 

obtain  the  zeros  of  g(x)  .  The  o-.'s  are  set  to  0  except  a     „  in  case 
^  J  ^   n+2 

the  zeros  of  f (x)  are  obtained  to  arbitrary  precision. 

Also,  there  are  occasions  when  the  zeros  of  g(x)  =  a  „  x 
'  ^  n+2 

+  f(x),  »    small  in  magnitude,  are  to  be  determined  wherein  somehow, 
the  zeros  of  f(x)  have  been  previously  known,  solved  exactly  or  esti- 
mated.  For  instance,  in  servomechanisms,  guided  missiles,  and  elec- 
tronic circuits,  it  is  usually  expedient  to  simplify  the  differential 
system  by  neglecting  higher-ordered  derivatives;  then  frequently, 
after  studying  the  truncated  system,  it  is  necessary  to  reintroduce 
the  higher-ordered  derivatives  one  at  a  time  and  study  the  augmented 
system.  Here  we  can  take  f(x)  to  be  the  characteristic  function  of 
the  truncated  system,  and  g(x)  that  of  the  augmented  system,  that  is, 
the  differential  equation  f(x)  vath  a  leading  tenn  adjoined  correspond- 
ing to  the  next  higher-ordered  derivative  previously  neglected.  Using 
the  iterative  procedure  we  derive  here,  we  can  get  the  zeros  of  g(x) 
from  the  known  zeros  of  f (x) . 


CHAPTER  II 


CONTROL  OF  TRANSUTION  ON 
PENULTIMATE  REMAINDERING 


Before  deriving  the  iterative  formulas  for  the  zeros  of 

n+1       .  1 
ff(x)    =  fY    x"*  +  S  (a  +CV,  )x   ,  let  us  first  study  a  method  for 
S\^J        "^+2   ^       Q_  ^  i  i     ' 

n+1    .  -, 
determining  the  zeros  r^'s  of  ^  a^x   ,  which  can  be  used  as  input 

to  the  iteration. 

One  of  the  most  convenient  methods  for  determining  the  zeros 
of  a  polynomial  is  Lin's  Penultimate  Remaindering  [1],  which  is  an 
iterative  scheme  of  producing  one  factor  of  a  polynomial  at  a  time. 
Like  most  other  methods  though,  it  is  a  contaminating  procedure  as 
has  been  explained  in  the  introduction.  With  the  zeros  obtained  from 
this  as  input,  the  iteration  which  shall  be  derived  in  the  succeeding 
chapters  can  serve  as  corrective  means  and  give  a  more  accurate  set 

of  zeros. 

Let  us  now  analyze  the  convergence  of  Penultimate  Remaindering, 
specifically  whether  or  not  it  can  be  controlled  by  translation  of 
variable.  We  shall  restrict  our  attention  to  an  nth-degree  polynomial 
f(x)  with  real  coefficients,  which  for  the  purpose  of  convenience 
stands  for  either  the  original  polynomial  whose  zeros  are  to  be  deter- 
mined or  for  one  of  its  depressed  forms,  and  the  factor  being  sought  be 
either  linear  or  quadratic. 


Our  problem  at  this  point  can  be  explicitly  stated  as  follows: 
Does  there  exist  a  real  number  h  for  which  Lin's  Penultimate  Remainder- 
ing shall  converge  for  f(x^^h)? 

n+1    .  ,   n 
Let  f(x)  =  E  a.x  ~  =  n  (x-v. )  and  the  divisor  d(x)  = 

(x-P  )(x-P  )  >riiere  P  =  y     for  k  =  1,2.  From  [1],  the  eigenvalues  of 
the  matrix  of  errors  in  the  iteration  are 

Pk  =  -^l^^°k'^  "  V2  Pk"^  "  •••  '=2^k)'  ^  =  ^'2'  "^^^^ 

G(x)  =  f(x)/d(x)  =  x"-^  -K  c^  -  x"-^  +  ...  c.x  +  c,. 

n-^  £  1 

Similarly,   if  the  divisor  is  linear,   say  d(x)    =   (x-P^ )   where 

^1  "  ^1'   "^1  "  ~^1   '^''l'     ■*■  ^n-1  ^1~     "^    •••   ^2^1^   vit.ere 

n         .    - 

G(x)    =  f(x)/d(x)   =  E  c.x^"-^,   c     =  1. 

2.     1  '     n 

Of  course,  these  eigenvalues  were  derived  on  the  assumption  that  at 

each  stage  of  iteration,  the  deviations  of  the  coefficients  of  the 

trial  divisor  from  those  of  thr  true  divisor  are  small  enough  to 

allow  terms  involving  said  deviations  of  order  higher  than  one  to  be 

neglected. 

To  obtain  the  corresponding  expressions  for  the  eigenvalues 

p  (h)  of  f(x+h),  we  use  the  fact  that  the  zeros  of  f(x+h)  are 

■y.-h,  j  =  1,  ...  n  and  that  p.  may  be  more  conveniently  written  as 

n 

k     -c^  c,  „  o  "        n 

^  ^         (-1)"-'ttv.       nv. 

3  J      3  0 

for  d(x)    quadratic,    and  k  =  1,2. 


'x)   linear,  p 

=  1 

n 

TT 

2 

(Vj-B^) 

n 

2^J 

"   (Y.-h-P  +h) 
3       J         ^ 

n 

=  1 

n 

n 

_  3_ 
n 

(-rV 

n   (v  -h) 
3       ^ 

IT 

3 

(vj-h) 

So     Pj^(h)    =  1  - =  1   -  ^ for  d(x)   quadratic; 


^  (y.-PJ 
2       J     1 
Similarly,   for  d(x)    linear,   p    (h)= +  1. 

TT    (y  -h) 
2        J 

Our  objective  then  is  to  find  a  real  value  h  for  which  |p,  (h)|  <  1 
for  k  =  1,2  or  for  k  =  1,  according  as  d(x)  is  quadratic  or  linear. 

Case  1.  There  exists  P,  equal  to  some  y,,    j  =1,  ...  or  n 
^ d 

This  implies  that  for  any  h,  except  possibly  h  =  3  if  real 
P,  (h)  =  1.   But  h  cannot  be  P  ,  for  if  it  were,  f(x+h)  would  have  0 
for  a  zero  and  the  Penultimate  Remaindering  cannot  be  applied  in  that 
situation.  Hence  there  is  no  value  of  h  for  #iich  the  iteration  would 
converge  in  this  case. 

To  avoid  such  an  exceptional  case,  however,  we  can  first  obtain 
all  the  multiple  zeros  of  f(x)  and  leave  a  depressed  polynomial  with 
simple  zeros.   One  way  to  do  this  is  to  apply  the  Euclidean  Algorithm 
for  determining  the  greatest  common  divisor  or  f(x)  and  its  first 
derivative.  To  illustrate,  if 

r      P-i 
f(x)  =  TT  (x-e^)  ^    P^  >  1, 

r      p^-1 
f  (x)  =  n  (x-e  )      H(x)  where  H  is  a  polynomial  of  degree  >  0, 


r      p.-l 
then  g(x)  =  (f,f')  =  TT  (x-e.)     whose  zeros  are  the  multiple  zeros 

of  f  (x) .  Here  (h,k)  stands  for  g.c.d.  of  h  and  k.  If  g  is  readily- 
factorable,  we  can  determine  the  multiple  zeros  of  f(x),  including 
their  respective  multiplicities.  Otherwise,  we  may  as  well  consider 
^(x)/g(x)  as  our  original  polynomial  instead  of  f(x),  for  the  zeros 
of  the  former  are  simple  and  precisely  the  same  as  the  zeros  of  the 
latter.  The  multiplicities  of  the  zeros  of  f(x)  obtained  by  the  iter- 
ation on  f(x)/g(x)  can  be  determined  later. 

Case  2.  All  the  zeros  of  f(x)  are  simple 

n 

Let     TT   (y  -3  )    =  A     +  B  i,   k  =  1,2,   vrtiere  for  each  k 
3       J     k  k         k 

A,  and  B  are  not  0  simultaneously.  Let  the  zeros  of  f (x)  be 

\'^2'   ^2t-l  =  \   ^  V'  ^^  ^2t  =  S  -  V'  "^^^^  t  =  2,  ...  s; 

V,  =  c  for  J  =  2s+l,  ...  n.  Here  s  is  an  integer  within  the  closed 

interval  [1,  n/2]   and  the  c's  and  d's  are  real  numbers  such  that  d  >  0 

A  slight  modification  will  be  necessary  in  case  d(x)  is  linear,  namely, 

that  the  complex  zeros  be  Vpx  and  y^     ^  ,  t  starting  from  1,  ...  s; 

and  the  real  zeros  aside  from  v,  =  P.,  be  y.^J  =  2s+2,  . . .  n,  where 

se[0,§]. 

Three  lemmas  are  worthy  of  note  at  this  point: 

n 
Lemma  1.  tt  (v.-h)(or  its  corresponding  expression  \ft\en   d(x)  is 
3   «J 
linear)  is  always  real  and  nonzero.   If  all  the  zeros  of  f(x)  are  com- 
plex, the  said  expression  is  positive. 
Proof. 

n   (v,-h)    =  TT   [(c   -h)    +  d  i][(c   -h)-d  i]   TT   (c   -h) 
3J  t^  tt  t-jj 

=    ^[{c.-h)'^    -H    d^]     TT     (c     -h), 


which  is  clearly  real,  and  nonzero  provided  h  ^  any  zero  of  f(x). 
In  the  above  the  empty  product  =1.  If  all  zeros  are  complex, 

n  2    2 

TT  (v.-h)  =  T7  [(c,-h)  +  d.  ],  which  is  positive. 

3   J      t    ■t       ^ 

Lemma  2.  When  P.  's  are  real,  tt  (y  -^,)    (or  its  corresponding 
expression  when  d(x)  is  linear)  are  real  and  nonzero.  For  two  succes- 
sive real  P,  's,  they  are  of  the  same  sign,  positive  v*ien  the  P,  's  are 
the  two  least  real  zeros. 

Proof.  The  first  statement  follows  from  the  proof  of  (1), 

replacing  h  by  3  .  To  show  the  other,  let  us  for  the  purpose  of  con- 

n 

venience  let  tt,  stand  for  tt  (v_,-P,  )  .  Note  that  the  sign  of  n  depends 
k  3   J  k  °  K       " 

only  on  the  second  product  n  (c.-0,).  For  P,  to  be  the  two  least 
real  zeros,  each  factor  in  that  product  is  positive,  so  tt  is  positive, 
for  each  k.  When  the  P.  's  are  not  necessarily  the  least  but  still 
successive,  there  would  be  as  many  negative  factors  in  tt  as  in  tt 
hence  they  would  be  of  the  same  sign. 

Lemma  3.   If  the  P  's  are  complex  conjugates,  '"^.'s  are  complex 
conjugates,  too. 

Proof.  To  see  this,  let  P  =  c  +d  i  =  ?  .  Expand  rr  as 

I  [(v=i)^(V^i)^]f^v^i)^^-VVi^2s".it^=j-^i)-V^  = 

^  [(VC^)^  ^4-4-   2d,(c^-c^)i]  ^n^[(c.-c^)-d^i]. 
On  the  other  hand, 

^2  =  i  [(V^l)^^V^l^i^f^V=l^"^-V^l)i]  2s'll^^'j"'l^'''l^  ^ 

I   [(VC^)^  ^d^  -  d^  +  2d^(c^-c^)i]  J^^  [c.-c^)+d^i].  Hence  n^  =  7F^. 


Let  us  now  go  back  to  our  main  problem^  that  of  showing  the 
existence  of  a  real  value  h  for  which  |p,  (h)|  <  1,  k  =  1,2. 

Theorem  1.  Let  f(x)  be  an  nth-degree  polynomial  with  real  coefficients 
and  simple  zeros  and  having  at  least  one  real  zero.  Then  there 
exists  a  positive  number  h  for  which  Penultimate  Remaindering, 
applied  to  f (x+h) ,  will  converge  to  a  linear  or  quadratic  factor 
according  as  n  is  odd  or  even. 

Proof.  If  n  is  odd,  we  know  there  is  at  least  one  real  zero. 

Let  P^  be  the  smallest  such.  Then, 

1  n 


n  (Y,-h) 
2   J 


where,  as  pointed  out  in  Lemma  2, 


A  =  n  (v,-e,)  >  0. 
2   J  1 


Since  tt  (y.-h)  must  be  positive  to  have  |p,  (h)  |  <  1,  let  it  be  C  >  0. 
2   J  ^ 

|p^(h)I    =    |l  -  ^1    <  1    <^  |C-A|<  C    <=>-  C  <  C-A  <  C    <=>C  >  A/2. 

n 

C  has  to  be  big  enough  so  that  the  equation  n   (y.-h)    =  C   would  have  a 

2   J 

real  solution  in  h.  For  practical  purposes,  we  can  estimate  h  provided 

we  know  an  estimate  P'  of  3,,  for  since  A  =  G(0^),  where  G  =  f(x)/d(x), 

A  can  be  estimated  by  g(P'),  where  g(x)  =  f(x)/(x-p')  -  G(x),  from 

which  C  can  be  estimated,  and  consequently,  h  can  be  estimated  by 

n-1 

V  C   provided  C  is  large  enough. 


10 


Henceforth,  let  us  assume  n  even.  The  rest  of  the  real  zeros, 
if  any,  may  be  obtained  in  a  similar  fashion  as  above,  either  singly 
or  in  successive  pairs  starting  from  the  smallest.  If  in  pairs, 

^(«  '  1  -  i =  ^  -  r 

Trjhere  A,,  A,„  >  0,  from  Lemma  2.  C  must  necessarily  be  positive,  too, 
to  have  |p,  (h)|  <  1.  Moreover, 

|pj^(h)|  <  1    '^=>   |c-Aj^|  <  C  <s=>-C  <  C-k^   <  c  <=$>C  >  k^/2,   k  =  1,2. 

<=^C  >  max  fAj/2}  . 

n-2 
So  by  choosing  C  large  enough,  we  can  estimate  h  by  V  C  . 

Thus  we  have  shown  that  a  real  zero  or  a  pair  of  real  zeros 
can  always  be  made  "stable"  by  a  proper  choice  of  h.   By  a  stable 
zero  is  meant  one  wherein  the  iteration  can  converge  to  its  corre- 
sponding factor.  Therefore,  all  the  real  zeros  of  f(x)  can  be 
obtained  using  Penultimate  Remaindering,  with  translation  on  x, 
if  necessary. 

Theorem  2.  Let  f(x)  be  an  nth-degree  polynomial  vdth  real  coefficients 
and  with  all  its  zeros  complex  and  simple.  Then  there  exists 
a  positive  number  h  for  which  Penultimate  Remaindering,  applied 
to  f (x+h) ,  will  converge  to  a  factor  corresponding  to  a  pair 
of  complex  zeros. 


11 


Proof. 

pj^(h)  =  1  -  -^  where  n  =  A  +  Bi  and  tt  =  A  -  Bi^ 
A,  B  real  and  C  >  0  using  Lemmas  1  and  3.  Since  p,  (h)  and  p„(h)  sire 
complex  conjugates,  |p, (h) |  =  |p  (h)|,  so  that  we  need  to  consider 
only  one  of  them. 

p^(h)  =  1  -  A^  =  (1  _  A/C)  -  (B/C)i. 

From  this,   it  is  clear  that  A  must  necessarily  be  positive  in  order 
to  have    |p^(h)|    <  1,   for  if  A  <  0, 

|p^(h)(^   =    (1-A/C)^  +   (B/C)^  >  1   since  A/C  <     0,    from  C  >  0. 

Using  Lemma  i;  below,    we  know  there  exist  two  zeros  of  f(x)   for  v^ich 

A  >  0;    hence  we   choose  these  two   zeros  as      P     and  0 

2  2 

|p^(h)I    <1    .^(l-A/C)^  -.    (B/C)^  <  1    <^1  -  ^  +  ^  +  2_  <  1 

«J=>C^   -   2AC   +  A^   +  B^  <  C^    <^C  >  Lll.   . 

ZA 

n-2 

Then  take  h  >  V  C   where  C  is  large  enough. 

Hence  a  pair  of  complex  zeros  of  f(x),  namely,  that  which 

corresponds  to  A  >  0,  can  always  be  made  stable  by  the  choice  of 

a  sufficiently  large  h.  By  -repeating  the  iterative  process  on  each 

depressed  polynomial,  we  can  thus  obtain  all  the  complex  zeros  of  f(x) 

Lemma  h.     If  f(z)  has  real  coefficients  with  all  of  its  zeros 
V.  complex,  then  there  exists  a  zero  &  vdiose  imaginary  part  is  posi- 
tive such  that  Re^  tt    (v.-P))  >  0. 
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Proof. 


^         (vj-P)  =    "      (P-V^)  since  the  degree  n  is  even, 


f(z 

) 

(z 

-3)( 

Z-'F) 

f( 

z)   - 

jf(e) 

z   - 

B 

1 

X 


z=3   ^ 


r'(g)   -2if'(p) 

2i  Im(P)    Im  (P) 
-2i[Re(f'(P)  +  iIm(f'(P))] 

uttfi 


So      Re(  n    (y  P))  =  ^^^I'liV 


Taking  P   such  that  Im(P)   >  0,  Re(     n     _  (y^-P))   >  0  if  and  only  if 


v>P,F  ^ 


Im(f'(P))  >  0.  To  show  this  is  true  for  some  B,    consider  the  function 

■jTT — r  ,  and  the  contour  integral  of  the  function  around  the  curve  C, 

composed  of  the  real  line  and  the  upper  semicircle.  C,  therefore, 

encloses  half  of  the  zeros  of  f(z),  namely,  those  whose  imaginary  parts 

are  positive.  By  Cauchy's  Residue  Theorem,  j    dz/f(z)  =  2tt  ±T.   r.  where 

C  -^ 

r,  stands  for  a  residue  of  l/f(x)  at  its  pole  inside  C,  and  the  summa- 
J 

tion  taken  over  all  such  poles.  Since  the  poles  of  l/f(z)  are  pre- 
cisely the  zeros  of  f(z)  i^ose  imaginary  parts  are  positive, 

r.  =  lim  (z-Y,)(l/f(z;)=  [  lim  iM]"^  =  ^^^  , 

if  y     is  such  a  zero.  Note  that  f ' (y  )  /  0  since  y     is  a  simple 
J  J  J 

zero  of  f (z) . 
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Therefore   J  dz/f(z)  =  2ni     E     l/r'(v  ) 
C  lm(v.)>0       J 

.  2TT  .V  [Re(f'(Y.i))    -  llm(f'(vj))] 
k'(v^)l^ 

/lm(f'(v,)\  /Re(f'(v.))i\ 

=  2E  i^  +  2E| J_      . 

lF(vJ?  |f'(Yjl' 


On  the  other  hand, 

J  dz/f(z)   =    lim    1/  dx/f(x)   +    f  ^l1^ 
C  R-»  L-R  0  f(Re^'f') 

The  first  integral  is  a  positive  value   since  as  x  -♦  +  <»,   f(x)  -  ", 
using  n  even,   \ih±le  the  second  is  0  since  the  denominator  of  the  inte- 
grand is  dominated  by  R     vrtiere  n  >  1.     Hence  J  dz/f(z)   >  0.     Comparing 

C 

this  with  the  first  result,  we  deduce  that  in:m(f'(v  ))  >  0,  from  which 

J 

we  further  deduce  that  there  is  one,  say  3,  of  those  v  ' s  such  that 
Im(f'(P))  >  0.  This  proves  the  lemma. 

Concluding  Remarks 

We  have  thus  proved  that  Lin's  Penultimate  Remaindering  can  be 
made  to  converge  to  a  quadratic  factor  by  a  proper  choice  of  some  posi- 
tive number  h.  It  would  be  very  nice  if  h  can  be  estimated  which  will 
work  for  all  the  zeros.  However,  this  is  very  difficult  to  determine, 
if  not  impossible.  For,  even  if  we  have  all  the  estimates  of  the  zeros 
we  are  seeking,  still  we  have  no  way  of  finding  out  if  a  certain  h  is 
good  enough  to  make  a  particular  pair  of  zeros  most  stable,  in  the 


Ik 


sense  that  the  said  zeros  would  be  the  first  ones  to  come  out  in  the 
iterative  process,  because  if  they  are  not,  then  further  translation 
might  still  be  necessary  on  the  depressed  polynomial. 

EXAMPLES 

Example  1.   f(x)  =x^-x^-x  -x-2=0 

For  h  =  0,  the  following  machine  output  is  given  for  the 
coefficients  of  the  divisor  x  +  b(2)x  +  b(l) : 

b(l)  b(2) 

2.0  1.0 

2.0  -3.0 

-0.6666667  -1.6666667 

-2.'^714287  -0.711;2858 

-I.I462686U  -1.2686565 

-2.ii890507  -0.755U76U 

-I.533367I4  -1.2333116 

-2.U357176  -0.7821669 

-1.5806082  -1.2096235 


It  is  not  converging. 

With  h  =  1.5,  the  iteration  converges  although  not  rapidly: 

f(x+h)  =  x^  +  5x^  +  8x^  +  2.75x  -  U.0625 

-0  5078125  0.3U37500 

-0.5881521  0.7Uoii567 

-0.7li75877  0.9670812 

-0.8380733  1.189232U 

-0.9U3U391  1.3803033 

-1.0292210  1.5618733 

-I.IIOI8U1  1.7185219 

-1.170lii;83  1.81190214 

-I.211I1O62  1.9222638 

-1.2328593  I.966OI6I4 

-1.2U31159  1. 98607 08 

-I.2I472219  1.99U5306 

-1.2i;895U2  1.9978617 

-1.2ii9578l  I.9991815 

-1.2U98h72  1.9996808 
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From  the  zeros     of  x^  +  1.9996808  -  1.2li98ii72   and  the 

depressed  polynomial,  we  obtain,    after  adding  1.$  to  each  of  them, 

the  zeros  of  f (x) : 

2.0  as  compared  to      2.0 

-0.99968  the  actual  zeros  <-1.0 

-0.00016  ±   0.9998Ui  +  i- 

Example  2.  f(x)  =  x^  -  2.075x^  +  .668x^  -  .091x  +  2.28 
For  h  =  0,  the  following  output  shows  divergence: 

b(l)  b(2) 

3.U131736  -0.1362275 

-0.75765143  -2.1687U13 

1.3996702  -0.0122631; 

-3.0120230  -3.6938910 

0.23602U1  0.1;9535ii3 

1.3370773  0. 3021^056 

U5. 7 2U91h6  61.92U6697 

0.0005819  0.7U68639 

0.8216316  -0. 032201a 


Taking  h  =  2,    so  that  f(x+h)    =  x^  +  5.925x-^  +  12.2l80001x 

+  9.68IX  +  Ii.l7,  the  iteration  converges: 

0.3ia2997  0.7923555 

0.5339ii31  1.0152895 

0.622^550  1.0537692 

0.6152719  1.0288571 

0.6380735  0.9979128 

0.6258326  0.98109U9 

0.6185358  O.9770U08 

0.6163983  0.9786261 

0.6167771  0.9809365 

0. 61761487  0.9822561 

0.6182072  0.982626U 

Thus,   from  the  zeros  of  x^  +   ,982626lx + .618072  and  the 

depressed  polynomial,   we  derive,    after  adding  2  to  each  of  them,   the 

zeros  of  the  original  polynomial,   namely, 

1.50869  ±  .61386  i 
-0.1i7119  1  .79782  i. 
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Check:   n(x-r*)  =  x^  -  2.075x^  +  .6679878x^  -  0.090i426lx  +  2.2776718. 

Remark:   Since  the  coefficients  of  TT(x-rT)  differ  from  those  of  the 
original  by  significant  amount,  we  can  improve  the  zeros 
obtained  here  by  applying  the  iteration  in  the  next  chapters, 
(See  Example  3  of  Chapter  IV.) 

Example  3.  f  (x)  =  x"^  -  8x^  +  lipc  -  12 

For  h  =  0,  trying  for  a  linear  divisor,  the  iteration  fails 

to  converge: 

b(l) 

-0.857lii29 
-1.5233160 
-2.9027832 
15.0731962 
-0.0331687 
-1.5153726 
-2.97U688i; 
12.6ii8U070 


Taking  h  =  5,   however,    a  linear  divisor  is  obtained: 


-1. 

-0.6591671 

-1.2100789 

-0.8978157 

-1.0565051+ 

-0.9707772 

-1.0156622 

-0.9917623 

-1.00ii3762 

-O.997687J4 

-1.0012255 

-0.9993515 

-1. 00031^  3U 

The  zeros  of  f(x)    are  then  found  to  be  6.0003ii  and   .999^3> 

+  1.00086i   (as  compared  to  the   actual  zeros  6.0  and  1  +  i) . 
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Example  h.      f(x)  =  x^  +  2^6. 8x^  +  l5l.9x^  +  258. 5x^  +  88. 5x^ 
+  k.Sx  +    .OS 

Applying  Penultimate  Remaindering  vdthout  translation,  two 
zeros,  -.01579  and  -.Oli28,  are  obtained  after  only  three  iterative 
stages.  However,  for  the  depressed  polynomial 
x^  +  256.7Ulii093x-^  +  136.8572083x^  +  250.3082123x  +  73.7ii23019, 
the  iteration  diverges. 

Applying  then  translation  of  variable  with  h  =  3,  the  poly- 

nornial  is  transformed  to  x^  +  268.7UlU055x^  +  2501. 529816 2x 

+  8111.U692993X  +  9069-3995361,  and  the  output  for  the  coefficients  of 

the  quadratic  divisor  sought  converges,  though  slowly: 

b(l)  b(2) 

3.6255U2  3.2ii2603l4 

5.51402665  U.  3670773 

6.760910U  1;.95U9237 

7.6359290  5.3278508 

8.3169955  5.5939965 

8.88lli500  5.8001015 

9.37314536  5-9698291 

9.8211602  6.1165862 

10.2i;39576  6.2i;86506 

10.656587h  6.371U687 

11.0711598  6.1;887165 

11.1980766  6.6026727 

II.9I457625  6.71UOI32 

12.Ul82l;27  6.8206988 

12.9077997  6.9152735 

13.3762U62  6.9789250 

The  zeros  of  x^  +  6.9789250  +  13.3762i;62  are  -3.1;89ii625 
±  1. 095141  from  iriiich  two   zeros  of  the  original  polynomial,   namely, 
-.Ii89ii625  +  1.095Ui  are  obtained. 
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2 
The  final  depressed  poljTiomial  is  x  +  26l.762li805x 

+  661.3328808  whose  zeros  are  -259.211lii02  and  -2.5513)402;  hence  two 
zeros  of  f(x)  are  -256.211lU02  and  .U1486598. 

Since  we  have  used  a  factor  whose  coefficients  have  not  con- 
verged yet,  it  is  possible  that  the  zeros  obtained  from  it  and  conse- 
quently the  succeeding  zeros  found  may  not  be  a  close  estimate  of  the 
actual  zeros.  Hence  the  iteration  process  which  shall  be  established 
in  the  n?xt  two  chapters  can  be  used  to  improve  those  zeros.   (See 
Example  10 of  Chapter  IV.) 

Example  5-   f(x)  =  x^  -  12x^  +  50x^  -  7Gx^  -  iOx^  +  2x  +  390 

As  in  previous  examples,  Penultimate  Remaindering  is  divergent 

for  h  =  0. 

For  h  =  6,  f(x)  becomes  x  +  2lpr  +  230x^  +  1130x-^  +  3019x^ 

+  i|0l46x  +  1950. 

The  following  output  shows  convergence,  but  quite  slowly: 

b(l)  b(2) 

0.6U59092  I.3UOI789 

1.1123276  1.9952l;ii9 

1.U510382  2.U011633 

1.7056U51  2.68303314 

I.90321UO  2.8925371 

2.0606095  3.055U308 

2.1887029  3.1861302 

2.29I47577  3.293U612 

2.3837912  3.3831297 

2.1i593760  3.I4590352 

2.52U1261  3.5239U91 

2.5800025  3.5799099 

2. 62850914  3.628I4607 

2. 67082141  3.6707982 

2.7078851  3.7078715 

2.7I4OI45I6  3.7l40i4i4l43 
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At  this  point,  instead  of  continuing  the  iteration,  it  might 
be  vd.se  to  consider  the  above  coefficients  as  those  of  a  divisor  of 
f(x+h)  -with  the  intention  of  applying  later  as  a  corrective  measure 
the  iterative  process  in  the  next  two  chapters,  since  the  latter  con- 
verges quite  fast  T'dien  the  input  r,  are  not  too  far  off  from  the 
actual  zeros. 

Applying  again  Penultimate  Remaindering  on  the  depressed  poly- 
nomial x^  +  20.259555x^  +  I5l.l798088x^  +  507.877880x  +  701^.1879913: 

U.6I487251  3.3527761 

7.8116007  ii.  7620618 

10.0788071;  5.5363917 

11.7581;51i7  6.002657U 

13.0062971  6.28U1839 

13.9031010  6.ii365376 

lU.h925766  6.ii977605 

li;.80li2977  6.U81i2725 

lU. 8710879  6. 1187196 

1U.7U7U96  6.2208961 

lii.U8i;8UU3  6.3232256 

lU. 1796360  6.1320639 

13.8976722  6.0698298 

13.6860651  6.0380195 

13.5619U7  6.0327021 

13.5182898  6.0U57873 

Again,   taking  x^  +  6.0l457873x  +  13.5182898   as  divisor,  the 

depressed  polynomial  is  x^  +  li4.213768Ux  +  52.0280985-       Using  the 

quadratic  formula  on  each  of  the  divisors,   the  following  zeros  of  f(x) 

are  determined: 

-1.0000222  +  6  =  U. 9999778 

-2.7U0l;221  +  6   =   3.2595779 

-3.0228936  +  2.099291  +  6  =  2.977106ii  +  2.0929i 

-7.10688U2  +  1.233i  +  6  =  -1.10688U2  +  1.2331- 

As  in  the  previous  example,  to  improve  these  zeros,  we  will  use 

them  as  input  to  the  iteration  which  shall  be  established  in  the  next 

chapters. 


CHAPTER  III 

ESTIMATES  OF  THE  ZEROS  OF  g(x) 

For  the  purpose  of  estimating  the  zeros  of  g(x) ,  we  propose  to 

derive  a  set  of  formulas  -vriiich  will  determine  the  values  of  Tl's  to  be 

added  to  the  known  r . ' s  to  obtain  the  zeros  r. 's.  These  values  can 
i  1 

also  serve  as  the  first  set  of  iterates  in  case  the  iterative  procedure, 
which  shall  be  discussed  next,  is  to  be  employed  to  get  better  esti- 
mates of  the  desired  zeros. 


""^         i-1 


Case  1.      t(x)   has  a  simple  zero 

Theorem  1.  Let  r, ,  i  =  1,  ...  n  be  the  n  zeros  of  t(x)  =  T.     a  x"^""^, 

and  r, ,  i  =  1,  . . .  n  and  r  ^  be  those  of 
i'      '  n+1 

n+1  .  , 

g(x)  =  cv^^2  ^    ^     I     ^h^  \U       . 

Writing  r,  =  r.  +11.  where  T],  —  o  as  the  o-  's  -  0,  for  i  =  1,  ...  n, 
^  1    1    i       1  q 

and  r "   =  -l/c  ^  +  "^     -, ,   then 
n+1      n+2    n+1' 

Q  K 

/i  \  Ti         Q=l  for  all  k  such  that  r,    is  a 


'k 


n   (r,  -r. 
L/k     k     X 


simple  zero. 


n  n 

(5)  \ 


n+1 


-11         n         n+2   ^      ^.        i  w         1     i'    ' 


w<i 
i,w=l 
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a         may  be  set  to  0  by  assuming  t(x)  and  f(x)  =  g(x)  -  ci         x 
to  be  monic  polynomials. 

Proof.  g(x)  =  a^^^   TT(x-r^)(x-rst^^)  =  ^  o  x*^"  +  "(x-r^) 

(6)  n+2   ^   i  i^'    a  ^     n+1' 

=  E  a  x^"-"-  +  rr(x-rj 
q  1' 

where  from  here  on,  i  ranges  from  1  to  n;  q,  from  1  to  n+2,  unless 

otherwise  specified. 

To  solve  for  \,   let  x  =  r  in  (6): 

(_T1,  )  n  {v.-r.A.){a     .  r, +1  -  a  ^  T\        )    ^  Ya  v}'"^    ^ 


^ 

TT  (r, -r.)[l l—][\  +  a     Ar,  A  ^^)] 

^j^       k  i'*-    ^k~^i  n+1  J 

Expanding  the  denominator  as  power  series  and  neglecting  terms 
in  Q"  '  s  and  11.  's  of  order  greater  than  1, 

-Eo  rj-^                          Tl                    Tl? 
11     =     ^    ,^^      ,     n      [1  +  -JL-  +  ^  +    ...    ]   X 

fi-V2^vVi)^  V2^vVi)'----  ^ 


Scv  r.^-^  ^ 


aJ^ assuming  l:.^^2^rj^-Tl^^^)  I    and 


(r,-r.) 
k     1 
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each  less  than  1,    which  are  valid  for  a   's   sufficiently  small,    since 
Tj.   -  0  with  Q-   's,    and  from  (6)    by  setting  x  =  0, 

(-1)   Q-T    +  rrr. 
^         ^  ^  '      1 1 

%+2     n+1       ~  ~     TT(r.    +  Tl.) 
1         1 

which  -»  0  with  »   '  s. 

q 

To  obtain   (5),      set  x  =  -^^^^2  ^"  ^^^*     "^^^^ 
\^2  "  °'n-.2  "  '^n+2  ""^  ^ 


Tl(l+r     c        )(-!)"                                           TT[l+(r  +T]   )cy        ] 
.  i-ill2 ^^       (_-ri        )(_!)"  ^     ^     "^2 

ff  -  n         n+2   n+1  »  ^  n 

n+2  n+2 

Multipljdng  by  a    (-1)   and  dropping  the  term  in  ff  ,  since 

Q-    =  0  from  f(x)  and  t(x)  monic, 

n 

u=l       n+l-u^  n+2'     ^     i  n+2'-' 

^S)     \+l  =  • 

n+2       n+2  i    i 

Expanding  the  denominator  as  power  series  and  dropping  powers  in 

0/   's  and  T|.  's  in  the  numerator  of  order  greater  than  2, 
q       1  b  J 

-1  2 

■"  [_i  -  a  a    +  1  +  Er.  a    +  E  E  r,r  c   ]  x 


n+1  -  a  ^       n  n+2        i  n+2   ,<,•   i  w  n+2 
n+2  w^i 

2    2 
n[l  -  a     ^(r.  +  T].  )  +  a-  _  r.  ] 
'■     n+2^  1    i'    n+2  i"' 


-1  2 

[-a  a     ^  +  Er.  (V  „  +  EEr,r  a-  „][l-Er.  o'  _] 

n  n+2     1  n+2    ^.   i  w  n+2-"-      i  n+2-' 


n+2    "  "  "     -  .^  -   ^^j_ 


~  [q-  -  Er.  -EEr.  r  a     _  +  (Er.)  o-  „] 
-   n     1    ^.   i  w  n+2     i'   n+2-' 

W<1 
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Case  2.   t(x)  has  a  multiple  zero 


Theorem 


_2.  Let  r  be  a  zero  of  multiplicity  n-,  >  1,  and  r  , 


s  =  n^+1,  ...  n,  be  the  remaining  zeros.  Let  r 

t  =  1,  ...  n.,  :  r  ,  s  =  n, +1,  ...  n:  and  r  ^ 

'  1'   s'      1   '      '      n+1 


be  the 


n+1  zeros  of  g(x) ,  where  r=r  +  T|,r  =r  +''1,  and 

\j  b    s    s    s 

r  ^  =  -^ —  +  T]    and  g(x)  is  as  previously  defined. 
n+1  ct     -    n+1     ^^    '  ^  ■' 

Then  a  monic  polynomial  Z  C         y     of  degree  n^  having  the  T]  '  s  as 

Q  k+1  i  t 

zeros  can  be  constructed  where  the  coefficients  C,  ^  are  determined  by 


(9)    S+1  =^  —K 


,       n+1  a     -,qlr^~     k  ,   r  ^ 


kl  TT  (r-r  ) 
n^+1 


,    (q-k)        ,  m        s'x=r 
q=k   ^^  '  m=l 


S-m+i^^-'"^': 


Here  TT'-'^-'(x-r  )    denotes  the  mth  derivative  of  (x-r  )  with  respect 
^   s  x=r  s 

to  X,  evaluated  at  x=r. 

Proof. 

A  polynomial  r  C    y  having  the  T]  '  s  as  zeros  must 
0  ^""^  ^ 

n^  -k 
have  coefficients  C  ,  satisfying  the  relation  C^^^^  =  (-1) 

multiplied  by  the  sum  of  the  products  of  the  ''l+'s  taken  (n,-k)  at 

a  time,  k  =  0,  ...  n-, -1  and  C   -,  =  1.  Indeed,  this  is  the  value  of 
'  1       "i  "*" 


2ii 


def.  .k   "^1 


(1/kJ)  n^^^(r)  where  TT^^^(r)  =  '  -^  [  n  (x-r-T]  )] 
t  t  dx*^  t=l     ^ 

evaluated  at  x=r. 

n+1  n    n 

g(x)  =  S  Q-  ,  x^  +  (x-r)     n  (x-r  )  =^ 
q=0  ^^^  s=n  +1    ^ 

(10)  g^^^(r)  =  E  a       qlr'^-^/iq-k): 

q=k  ^^ 

But  from   g(x)  =  c^^^  ^^"Vl)  n(x-r^-H. ) 

=  a  „(x  +  -  Tl  .)  tt(x)  ■n(x),  where  from 

n+2^    a   ^„    n+1'  +    c, 
n+2        ^         ^ 

here  on,  t  ranges  from  1  to  m ;  s,  from  n-+l  to  n,  and  as  in  the 

previous  theorem,  i  from  1  to  n,  and  q  from  1  to  n+2,  unless  otherwise 

def.  ^g^ 

specified;  tt(x)   =  n(x-r-Tl  )  and  tt(x)   =  '  (x-r  -T\   ) , 
t  ^       s  ^  ^ 

we  can  proceed  by  induction  to  show  that 

/■i-i\  [k]/    X  ,      V     /k-lv  „[k-l-m].    s   „[m],    .     , 

(11)  g^  ^(x)  =  o/^.^k     Z  (    )  TT^     J(x)  n^   (x)  + 

{a     ^x  +  1  -  a  ^1\        )     Z  (^)  n^^""^(x)  n^'"^(x). 
^  n+2        n+2  n+1'  ^=o  "  *        ^   ^  ' 

It  may  easily  be  verified  for  k  =  0,1,  so  assuming  the  identity 

true  for  k,  we  can  show  it  tme  for  k  +  1  as  follows: 

k=l 

r  [k]/  ntI      ,   k-  /k-1^  r„[k-m]  /  X  „[m]/  s  ,  „[k-m-l]/  s 

[g   (x)]   =  a       k     E  (    )[n^    -"(x)  n^  ^(x)  +  n       (x)  x 

n+^  m=0  "^   t        s  t 

„[m+l]/  \i   r       T       Tl   1  V  ^kv  r„[k+l-m]  /  >, 
s  ""^2        n+2  n+1  jn=o  ™  t 

^[™]^  \   „[k-m]/  N  „[m+l],  v,        ^  /kx„[k-m]/  s  [ml/  n 
rr^  -"(x)  +Y     (x)  tt'    '(x)]  +  a  ^JJT  (x)^  ^(x) 

s        t        s  '^■'■'^  m=0  ffi  t       s 
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m=0  "'  t 


^J^.l^  Jk-m.l](^)  Jm]^^) 


m 


=0  ^   t  s 


1^+2  m=0  "^  t        s 

using  the  lemma  below  on  the  first  two  groups  of  summation  in  the  deriva- 
tive, with  V  =  k  and  k  +  1,  respectively. 

V  /V-l>,r„[v-m]  /  s   [ml.  ^   ^[v-l-m],  \„[m+l].  ^  _ 
Lemma.   E  (    )[n'-    -"(x)  n"-  -' (x)  +  n"-     ■'(x)tt'-    ■'(x)  ■ 

m=0      t        s        t         s 

E  (  )  tt'-  ""'■'(x)  n*-  ■'(x)  for  any  positive  integer  v. 
m=0    t        s 

Proof.  Left-hand  side  = 

''i\^-l)nt^-'^^x)n^"^](x)  .  ''i\^-^)n[^-l-"](x)nt"^"l](x)  = 
m=0  "  t       s        m=0  ■"  t         s 

^  m  .  1   m-x  X. 


v-1 
S[( _ 

m=l     '"        '"   -^     ^  g  "    t  3 


'A\'-)^{l:l)]-^''-''^^^''h-)    -  Cl'-)-^''h:<)-^''h:^    - 


/V-lv    [0],   N„[V],   V 

~     t  s 


^;\v.nt^-](x)n[-](x)  .  (:;)nf^^(x)n[Ol(x)  .  0-^''^^)^^''^-) 

\Wt  s 

m=l 

right-hand  side.   (QED) 


T  m  t       s  t      s 

m=l 
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Now,  equating  (10)  and  (ll)  evaluated  at  r, 

n+1  (y      Q '  r  k-1 

^   Q+1        -    1   V  /'k-ls„[k-l-m]  ,  \„[m]/  v 

,    (q-k)       n+2     ^'  m  '^.       ^    '  ^    ' 

q=k   ^^  '  m=0     t         s 

^V2^^^-V2  Vl^   Vm^"    ^^)"  ^^)- 
m=0   t       s 

[k] 
Solving  for  ti"-  ■'(r)(m=0  of  the  last  summation), 


(12)  J^^(r)  =f{"?  <.^  .qlr^-V(q-k)l  -  <.  „  k  ''e\^-1)  ^^"1-]  ( 
t        L  q=k  ^''-'-  """^    m=0  "^  t 


r)  X 


nt-^(:)/[l..   (r-n   )]  -  E  (;^)n[k-]  (r)n[-l  (r)  /n(r; 
s     ;     n+<^    n+l     m=l  i"  t       s     J-'  s 


Expanding 

l/n(r)  =  l/n(r-r  )[1-T1  /(r-r  )]  and  l/[  1+a^  -(r-\  _ )  ] 
s  s     s     s  n+^    n+l 

as  power  series  and  neglecting  terms  in  o'  's,  T] .  '  s  and  of  C  ,  's 
■^      "  q  '  i  v+1 

n   (r) 
(i.e.,  1 ' s)  of  order  greater  than  1,  we  arrive  finally  at 


n+l  a   ^  nir^"^    k  ,   r  , 
E  .Sll_.  E  (l^)n[^](x-r  )    C^   ,(k-m), 
[k],  ,   q=k   (q-k)  J     m=l^        s^^=r  k-m+1^ 

n"-  '(r) . 

t 


n(r  -  r 


From  this,  (9)  follows  immediately. 


CHAPTER  IV 


A  NONCONTAKLNATING  ITERATIVE  PROCESS  FOR 
DETERMINING  THE  ZEROS  OF  g(x) 


Case  1.   t(x)  has  simple  zeros 

Theorem  1.  Let  t(x)  have  zeros  r. ,  i  =  1,  . . .  n.  If  r.'s  are  the 

corresponding  zeros  of  g(x),  that  is,  r.  =  r.  +  T). ,  i  =  1,  . . .  r 

such  that  fl,  -*  0  with  the  o-  's.  and  r  ^  =  -1/a  o  +  ^  tj  'then 
i  q  '  n+1      n+2    n+1' 

(A)  i}. 

r,    +  U    '^     -*  r.    as  J  "*  °°  for  all  k  such  that  r,    is  simple  zero 

(J)     * 

'^     -'  -^  where 


of  t(x),  B:^d  -1/^,2  '  \^1        "^  Vr 
(13)     T]  ^^^^'    =  ^  ^ 

lA  _^  J  >  0 

(iw      Vi*'  ' m 

V2'ni:i*v2(''i*V"" 

and  T),  '^  and  11  i      are  the  estimates  (U)  and  (5),  respectively,  of 
K        n+1 

the  previous  chapter. 

Proof.  Comparing  (13)  and  (lli)  with  (7)  and  (8),  respectively, 
of  the  previous  chapter,  we  see  clearly  how  the  iteration  formulas  are 
arrived  at,  namely,  by  replacing  all  the  T] '  s  in  the  right-hand  side  of 
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(8)  and  (7)  by  Tl^-^^'s  and  those  at  the  left-hand  side  by  T)^^'*'-'-^  '  s. 
Let  us  defer  the  proof  of  the  convergence  of  these  iterative  proce- 
dures for  sufficiently  small  a   's,  and  go  to  Case  2. 

Case  2.   t(x)  has  a  multiple  zero 

Theorem  2.  Let  r  be  a  zero  of  multiplicity  n,  >  1  of  t(x) .  Let  r, , 

t  =  Ij  ...  n,  be  the  corresponding  zeros  of  g(x),  i.e., 

r.  =  r  +  T),  where  T]  -•  o  with  the  a   's,  and  r  ,  s  =  n, +1,  ...  n 
t       t       t  q  '      s       1  ' 

be  the  remaining  zeros  of  t(x) .  Then  r  +  T)^^^  "*  r,  as  j  -•  "^ 

(1)  "■'■    (i)  k 

•where  T|  "^  's  are  the  zeros  of  the  raonic  polynomial   E  C,  ,   y  , 

the  coefficients  C  ^"^^   being  determined  by 

/  .   ,  V  n+1  a  ^,qlr^  k-1  ,    ,  /  ^  -,  \ 

(1^)   ^vJ^      ^    =       ^  Uv\^        -  c     ,k  E  (^-l)(k-l-m):   C,      ^J^l)   X 

k+1  ]q=k        (q-k)l  n+2     q  '  m   ^  ^        k-m 

(J  >0) 

n[™^(x-r  -Tl^-^^)  -   [1+cy     JrA     ,^^^)]   x 

^        s     s     ^x=r       "•       n+2^        n+1       ^  ■" 

I  (^)(k-m)l   C^,      ^^'^K^'^^x-T  -Tl   (J))  / 

ni=l  m'  k+l-m  s     s        ^x=r    / 

kl  TT(r-r^-Tl^^«^))[l+a^^2(^-Vl^^^)^  "^^''^  ^^ 
stated  before,  i  ranges  from  1  to  n  and  s  from  n, +1  to  n,  while 
C,_^^   's  are  the  estimates  (9)  obtained  in  the  previous  chapter. 

Proof.  As  was  already  noted  in  the  previous  chapter,  C,  ^  is 

[kl 
precisely  n'-  ■^(r)/k!;  comparing  (12)  with  (1$)    above,  the  derivation 

of  the  iteration  formula  is  clear. 
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We  shall  now  establish  the  convergence  of  (13),  (ll^),  and 
(15)  simultaneously  for  sufficiently  small  values  a   's.  Without  loss 
of  generality,  let  us,  for  the  sake  of  simplicity,  assume  that  t(x)  has 
one  multiple  zero  r  to  which  corresponds  r  's. 

We  shall  presently  aim  to  obtain  a  matrix  equation 

e(J+1)  =  ^   E^^'^-'-^  +  N  E^'^^  i^ere  M  and  N  are  are  (n+1)  by  (n+1) 
matrices  to  be  determined  and  E  '^'^  is  a  column  vector  with  elements 


\ 


(J)  ^  ^^i^)   _  Cj^  for  k  =  1,  ...  n^  while  E^^^^^  =  T\^^^^   -  T\^ 

for  k=n-.  +1,  ...n  +  1. 

Let  us  then  go  back  and  consider  Case  1.  From  (7)  and  (13), 
for  k  =  n,  +1,  ...  n, 

(16)  E,^^""^^    =   -     E     c^  rj"-^]  n   (r^-r. -TIJ  [1+a     Jr-T]     J]    - 
k  q=i     q  k      ]ijiy^     k     i     i^*-       n+2^   k     n+1^  •' 

TT   (r,-r.-Tl.^'5))[i+^        (        Tl  E     ,^^h]     / 

j_^jj     k     1     1  n+2     k     n+1         n+1  |/ 

TT  (r,-r,-11,  ^^^)  n   ( r  -r, -Tl . )  [  1+a     „(r, -11     J]  x 
j_^jj     k     i     i         i/k     k     i     1  n+2     k     n+1 

[1+cy     „(r, -Tl     ^    -  E     ^^^h] 
*■       n+2^   k     n+1         n+1       ^ -* 

=  -     La  r,*^"^  <*   n   (r, -r.-TlJ[l+cy     ^(r, -T]        )]- 
^  ^       liA  n+1  ■■ 

n^-1 

,5,^V^s-\-  ^s^^^)t   ^^^k-)'  ^+1  ^  ^0     ^^k-)\+l^^^)^ 

f^%+2(vVi  -  \+i^^^)|/ijk^w\)  - 
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n^-1 


(j) 


^•^  n+2^^k  n+1^-"-^  n+2^^k  n+1    n+1  '^' 
where  from  now  on,  v  ranges  from  0  to  n, ,  unless  otherwise  specified. 
Neglecting  terms  in  E  "^  '  s  of  order  >  1,   numerator 

N  ~  -Sa  r^-^  I    n  (r, -r. -11,  )[l+o  Jr.-T\     J]  - 
~    '^  ^       I  i/k  ^     i  i     ""^2  k  n+1 


[  n  (r, -r  -T)  )  -  E  E^^^  tt   (r  -r  -Tl  ))]  x 
s/k  k  s  s^   p^k  P   s/p,k  k  s  s  ■■ 

n^-1 

[E(r, -r)"^  C  ,  +  E  (r, -r)""  E  ,^^^]   x 
■■  ^  k  '^   v+1    0    «^     v+1   ■• 


[l+(>  ^(r, -Tl  J  -  a  ^  E  ,^^^]'> 
'•   n+2^  k  n+1^    n+2  n+1   ■' 


where  from  here  on,  p  ranges  from  n, +1  to  n  unless  otherwise  specified, 
and  empty  product  is  defined  1. 

Expanding  and  neglecting  further  terms  in  E  "^  '  s  of  order  >  1, 
the  terms  independent  of  E  '^'^'s  cancel  out,  leaving 

(17)  N  ~  Ecy  r'^-^i  [E(r  -r)""  C  ^]{lw     Jr, -f\     J]E  E^^^x 
q  k   j  •■  ^  k  '   v+l"'^   n+2'  k  n+1'  p^^  P 

(  TT   (r  -r  -Tl  ))  -  n  (r  -r  -T]  )[l+a   ^r  -Tl  J]  x 
s/p  k  '^  ^  ^     ST^k  if  s  s     n+2  k  n+1 

n^-1  I 

E  (r,-r)^E  ,^^^    +  a     ^     n  (r  -r  -Tl  )  [E(r  -r)\  JE  ,  ^^H 
0    k  "^   v+1      n+2  s^k     ^     s  s'""  ^  k  '  v+l""  n+1   f 
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Similarly,  denominator 

n^-1 


2/       -n  NT.  -    /    T,    V.2 


D  ~  TT  (r  -r  -Tl  )ri+Q'   (r  -Tl   11 


n    k  '  v+1 
[1.^ ]x 


t" 


k  ^  v+1 


n   pA  ^    s/p,k ,.     n+2  n+1      ^ 

sA  '^  s  s  n+2^  k  n+1  j 

Expanding  the  reciprocal  of  each  factor  within  the  brace  as 
power  series,  and  noting  that  the  numerator  has  no  term  independent  of 

E  "^  '  s,  we  only  need  the  constant  of  the  resulting  expansion,  namely, 
1,  for  orders  of  E^'^'^'s  >  1  are  being  ignored.  Thus 

Varying  k  from  n, +1  to  n,  we  thus  derive  from  (17)  and  (18)  the 
(n-.+l)st  to  nth  row  elements  of  our  desired  matrices  M  and  N.  Note  that 
the  (n^+l)st  to  nth  rows  of  M  are  filled  with  O's. 

Moreover,  from  (8)  and  (lU), 

n 

-1  +  E  (-a     J^  a         ,  +  nd+r.Q-   ) 
(j.l)        1    "^^    "-^^^       ^  "-^ 

'"f^%.2^V\^'^^  -n[l.c.^^2^r.+Tl^)] 


n[l+a     Jr.+T][^hMl+o  ^^(r.+T].)] 
n+2     1     i  n+2     i     x 
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Numerator  N  vjithin  the  brace 


=  rr[l+a     ^(r  +T1   )^a     ^E   ^^^][E(-a     .)   "^      (l+o     -r)^C     , 
n+2     s     s       n+2  s       ^  n+2  n+2         v+1 

n, -1 
1  n.,  -V  /  .V 

0     n+2^    ^   n+2  ^  v+1   -'    '^   n+2^  i  i^-" 
Neglecting  terms  in  E  "^  '  s  of  order  >  1,  and  combining  terms, 


we  get 


n^  -V  ,    . 

-     n+2'-  n+2  n+2    '      v+l""  p         st^d         ^^"^2      s     s^ 


n^-1 


n,  -V  ,  .  ^ 

1     M.^       >,^v^       (j) 


+  tt[1+q'       (r  +11   )]   E     {-a     .)    ^      (1+a     .r)  E     , 
"^        n+2     s     s   -^    0  '^"*"2  n+2         v+1 

Similarly,  denominator  D  within  the  brace 

2 

~  TT   [l+cy     ^(r.+Tl.)].     Hence 
-  n+2     1     1 

n 

-l+E(-a        )^a         ,+n(l+r.Q'     „)rn^-l 
/..nx  .        n+2'     n-u+1  1  n+2'rl  n-,-v 

(19)E^^"^U— ^ ^ ^     ^-^n+2)  (l%+2^)''^ 

""^  t>     V[l+a     Jr.+Tl.)]  p  "+2  n+2 

n+2      '■        n+2^   i     i'-"  [ 

from  which  we  get  the  last  row  elements  of  M  and  N.     Note  that   (19) 
does  not  also  give  any  nonzero  contribution  to  the  matrix  M.     Note 
further  that  each  term  of  the  numerator  in  the  expanded  form  of  (19) 
has  at  least  an  order  -2  factor  of  o"   's,    so  that  each  coefficient  of 

q 

E  .  "^      and  E  '^     which  is  an  element  of  the  matrix  N  goes  to  0  with 


P 


the  <y   '  s . 

q 
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Leaving  this  for  the  meantime,  let  us  obtain  for  Case  2  an 
analogous  expression  for  E  '^  shen  k  =  1,  .  .  ,  n,  .  From  (l5)  and 
its  corresponding  expression  for  C  ^,   for  k  =  Ojl,  ...  n, -1,  / 

f"*l      a   k 

n(r-r-Tl^(^'][l«„,j{r-Vi'^'>lU:''('-V^l^'>== 


Proceeding  in  the  same  manner  as  in  the  previous  case,  that  is, 
expanding  out  and  neglecting  powers  of  E^^^*^  higher  than  1,  the  terms 
independent  of  E  ''•^^'s  drop  out,  leaving  finally  a  linear  combination 
of  the  E/''^^  's,  namely. 


(20)  E^^^(J^l)  ,     -  ^ir-r^-\)[l^a^^^ir-X^^)]a^Jf&-^^U 


(k-l-m):   rrf'^^rx-r  -Tl   )  E        ^^^^^ 

s     s'x=r     k-m 


3h 


,2  ^   ,k> 


n(r-r^-n^)[l..^^p(r-^^^^)]^S(-)(k-m) 
J"^J(x-r  -T]    )         E,    ,      (J-^1)    .n 


s     s'x=r     k+l-m 


(r-r  -T]   )  X 
s     s 


k-1 


k-1, 


^^%.2(-Vl)]V2  ^     g  (';;')  (k-l-m)!    C^_^  X 


[i:Ep(J)(^n^[-^(x-r^-'OJ_]   .n(r-r^_^Jx 


s     s  x=r- 


|2  ^   ,k 


(j) 


f^%.2^-Vl)]^(.)('<-)^S.l.J^E;^^x 


[m] 


^sV    ^^^-V\)x=r]^f[^-n.2(-\.l)] 


n+1  jf_]^ 


^1.  „  TTl^"^(x-r  -T|   )        ]    _   [i+cy        (r-T\        )] 
k-m  s     s^x=r-'        '•       n+2^       n+1'^-' 


[m], 


SnCk-m)'    n  nL™J(x-r  -Tl    )        }   EE 

■j_  tn  k+i-ra  s     s  x=r 


(j) 


n+1 


^3?p^^-V\»^  -n.2  "(-^s-\)[   ^  Vl"'"'^'/ 


q=k 


3^ 


(q-k)i -.^^2kY(';')(i^-i-)iv^>^(x 


-r  -Tl  )   ]  E  ^,^^H/ki  "^(r-r  -11  )  x 
s  s  x=r-'  n+1   (/         s  s^ 


n+2    n+1 

Varying  k  from  0  to  n, -1  in  (20)  we  can  get  the  first  n, 
row  elements  of  the  matrices  M  and  N.  This  completes  our  determina- 
tion of  the  matrices  M  and  N. 

Note  that  in  the  vector  equation  E^*^"^"'-^  =  M  E^'^'^'^^  +  N  E^'^^ 
the  nonzero  elements  of  M  are  found  only  in  the  first  n^  rows  to  the 
left  of  an  excluding  the  diagonal.  As  for  N,  note  that  each  of  its 
elements  has  in  every  term  either  a  factor  a  or  C,  of  power  greater 
than  0  and  since  C,  -  0  with  the  a  's,  together  with  the  fact  that 
the  T] '  s  ~*  finite  number,  we  can  make  each  of  its  elements  arbitrarily 
small  with  the  a   's. 

q 

To  illustrate  it  for  n  =  Ii  and  n,  =  3,  see  the  following  pages. 
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M=- 


2a, 


[r-r^-l|^J[l«,(r-ll5)) 


2a. 


Il-6''-"5*J   *  '^^^V 


N  = 


6       n  n    3  p 

-Sa  r?"V"(r, -r.-10.)    [1+ 
]_  q  1;     'i     U     1     1 


".(^Tlj)] 


6  ,3 


'T'-'  <""6 


-\'  tl*«6 


.[-l+E(-a   )^cy       +n(l+r,x 
"^       ]_       6       5-u  i 


2/_2 


Oi^\ay^\\^^{vf\^\ 


[-l+EC-ffJ^a^     +n(i+r.   x 
1       6       5-u  1 

2 

Q'^)]o'^(l+Q'^r)/n   [l+o-^ 


-[-l+E(-(y^)' 


(r.+Tl.)] 
1     1 


^\)^ 
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(K.i^')/^-vV'f^ 


^^(r-n^)]' 


^6(k.l^')/f^-V\^  ^ 


[l«v^(r-ll^)]' 


,q-i 


[r-r^-Tl^.l]/[r-r^-Tl^]2  x 
[l^^(r-Tl^)] 


(r-Tl^)]^  j 


[l-H>^(r-Tl^)]/2(r-r^-Tl^)2 
[l+c.^(r-Tl^)] 


(C2(r-r^-Tl^)-2C^)]/2(r- 


^r^/j^v^i 


(V^,)] 


o^     +n(l+r.x     I    [1-+E(-Q',)^Q'^     +n(l+r.Q'   )>i 
5-u  1       I  1       u       i>-u  1  c 

TT^[l+c.   (r.+Tl   )] 
oil 
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e(J-1)  =  M  E^J-1)  ,  N  E^J)  ^  E^J'l^  =  (I-M)-^^J\ 

Note  that  I-M  is  a  lower  triangular  matrix  with  diagonal  elements 
1  since  M  is  subtriangular  and  hence  (I-M)   is  also  a  triangular 
matrix  with  diagonal  elements  1.  Moreover,  each  element  of  (I-M)~'TJ 
consists  of  finite  sum  of  terms  each  of  which  has  for  a  factor  an 
element  of  N.  Therefore  each  element  of  (I-M)  "TJ  can  be  made  arbi- 
trarily small  in  magnitude  with  the  \(y   \'s.     Thus  for  sufficiently- 


small  magnitudes  of  »  ,  max  [r|p  |]  <  1  v±iere  (I-M)'"!]  =  ( 


P,.„,'' 


Therefore  the |eigenvalues|  of  (p  )  are  all  less  than  1.  This  proves 
the  convergence  of  the  iteration  scheme  established. 
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ILLUSTRATIVE  EXAMPLES 


A.  As  a  Method  for  Determining  the  Zeros  of  a  PolynoTuial 
from  the  Known  Zeros  of  Its  Lower-Degreed  Part 


Example  1 

Suppose  the  zeros  of  h(x)    =  x^  +  99.999992lx     -  800.ii9995U2x 
+  lU00.9990082x  -  1199.79li521  are  to  be  determined  vdierein  99 •  999992Ux'* 
-800.I;9995U2x^  +  lli00.9990082x  -  1199.79U521  had  been  previously 
studied  and  its  zeros  approximated  to  be  6  and  1  +  i.     We  may  consider 
instead 

where 

f (x)  =  x^  -  8.005001x^  +  1U.0099911X  -  11.9979l;65, 

and  apply  the  iteration  on  g(x). 

Here 


a^   =      .0020536 

a^   =  .0099909 

a     =  -.0050001 

%=   .0 

a^  =     .01 

I  by  the  machine. 

r  's:  1  +  i 

6.0 

\ 

\ 

-.000^385  +  .00268781  -.1^939229  -7.5050001 

-.OOC560U  +  .0027260i  -.U3li9750  -7-5683575 

-.0005656  +  .00269231  -.h3U72Ul  -7. 56141376 

-.0005655  +  .00269231  -.U3U7l;39  -7.56U1173 


ho 


Therefore,  the  zeros  of  g(x),  and  consequently  of  h(x),  are 

.999I43U5  +  1.00269231 

5.5652561 
-107.56U117U 

as  compared  with  the   actual  zeros  of  h(x): 

.999U3l;5  +  1.00269231 
5.5652561" 
-107.56iill7ii. 

Example   2 

Let  h(x)    =    .OOlx^  +  x^  -  iix^  +  8x^  -  8x  +  U 

Suppose  the   zeros  of  x     -  Ipc     +  8x     -  8x       It  are  previously 

known  to  be  1+i,   each  one  of  multiplicity  2.     To  determine  the   zeros 

of  h(x): 


r.'s: 

1+i 

1+i 

Q".  's: 

0. 
0. 

0. 

0. 

0. 

0.001 

CO^p 

co,2 

.OOli 

+  .OOli 

^°ir   '^°12  ^°21'   ^°22  \ 

.001    +  .0011        .OOli     +  .OOli       -3.9920000 
.000995  +  .00099651    .0039801  ±  .00099751   -3.9919533 

Hence,  the  iteration  converges  after  only  two  iterations,  giving 

the  zeros 

r!"'s:         -1003 .  9919U3ii 
^  1.0123783  +  0.96l48925i 

.9836U15  +  1.03Uli 

as  compared  to  the  computed  zeros   [}], 

-1003.9920iiO 
1.0123775  +  .96ii89lil9i 
.9836I425  +  1. 031422991. 


i;l 


B.  As  Corrective  Measure  for  Results  Obtained 


by 

Other  Methods 

Example  3 

f(x)  =  x^  -  2.075X 

3. 

.668x^  -  .091x  +2.28 

(from  Example  2  of  Chapter 

II) 

r.  '  s: 

1.50869  +  .613861 
-.U7119  1  .797821 

<y .  '  st 

.0023281 

.0005739 

.0000122 

.0 

.0 

.0 

1'  '2  '3'  h 

.0001738  +  .00018811         -.0001738  +  .00029i;Oi 
.0001733  +  .00018831        -.0001733  +  .00029391 

After  two  iterations,  the  zeros  are  found  to  be 

1.5088637  +  .6lU0h83i 
-.U713638  +  .79811391- 

Check:  n(x-rj)  =  x^  -  2.075x3  +  .6679998x^  -  .0909998x  +  2.2800001. 

Example  U 

f(x)    =  x^  +  200x3   -  I600x^  +  28OOX  -   2UOO 
(from  Flanagan-Maxfield' s  example  [3]) 


r^'s: 

.999615  +  1.001931 
5.75077 
-208.0 

Q-.  's: 

1 

-3.97903U1 
3.I4766235 
-1.5002289 
-  .25 

0. 

0. 

1|2 


'l'  '2 


-.0000096  +  .0000051   .0161016 
-.0000096  +  .0000051   .0161192 


.2339176 

.2339000 


Therefore,  the  zeros  are: 


.999605U  +  1.0019251 
5.7668892 
-207.7660999 


as  compared  with  actual  zeros,  [3] 


.993^^"^^   +  1.0019251 
5.766889  " 
■207.7660998. 


Example  5 

f(x) 
(from  Lin's  example  In  his  19ill  Iterative  process  [7]) 


x^  +  10.65x^  +  129x^  +  203. 5x  +  70 


estimates  r. ' s: 

1 


-1.3775 
-0.3775 
-ii.U55  +  9.651 


\ 

.0996209 
.1109277 
.1123392 
.1125191 
.1125119 


cy.  's  are:    11.2551U32 
^  O.6O523I1I 

-0.1265793 
-0.015 
0.0 
0.0 

\         h-  \ 

-.1003151  .00781i71  ~  .0070U05i 

-.1116175  .0078i;89  +  .007051i8i 

-.1130373  .OO78U90  +  .00705571 

-.1132171  .0078ii91  +   .00705581 

-.II32UOO  .OO78U9I  +   .00705581 


Therefore,   the  zeros  r. ,      are: 
'  1'  s 

-1.26i;958l 
-0.h907UOO 
-ii.UI;7l509  +  9.6ii29iii|ll. 

Check:     n(x-rj')    =  x^  +  10.6ii99999x^  +  128. 99999l43x^  +   203.ii9997lipc 


+  69.9997603. 
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Example  6 

f(x)   =  x^  -  6x^  +  U7x^  -  I8x  +  290 
(Lyon's  example  using  Lin's  19i|l  method   [7,10]} 


estimates  r.  ' s: 

-0.2855  +  2.5ii5i 
+3.2855  +  5.78i 

the  a,  '  s   are: 

.0936661 
-.m37273 
-.009ii032 
- . 0000001 
0. 

11  ,  Tl 
1'   2 


-.0010757  ~  .00211971 
-.0010762  ~   .00212071 


Therefore  the  zeros  are 


W 


.0010758  +  .001ii573i 
,0010762  +    .OOlU580i 


Check:      n(x-r'! 


3.28iiJ4238  +  5.77787921 
-0.28UU237  +  2.5h6U579i 

x^  -  6.0000002x^  +  i;6 . 9999995x^  -  l8.0000036x 

+  289.9999809. 


Example  7 


f(x)   =  x^  +  25^.79x^  +  152. Olx^  +  258.ii3x^  +  88.37x^ 
+  J4.52X  +   .0Uli3     (Luke's  example  [9]) 


r.  '  s: 

1 


a .  ' s  are : 

1 


-.OU59 
-.0130 

-.107  +   .9i;65i 
-.316  ~ 
-256.2 

;  .OOOU693 

.001^2755 

.10ii3ii82 

.1561;789 

.1270771 

.  0011021; 
0. 
0. 


v\ 


Check:  TrCx-r?')  =  x°  +  256.7899933 


x^  +  152. 


kh 


-.0006630    .0000385   -.001060  + 
-.0006629    .0000378    -.001061  + 

.00013771 
.00013771 

.0003Ui;3 
.  0003ii50 

- . 0006102 
-.0006102 

Therefore  the  zeros  r.'s  are: 

1 

-.U65629 
-.0129622 
-.1071061  + 
-.3156550 
-256.2006035 

.9U66377i 

x^  +  258.1;29977Ux' 


+  88.36999ii2x    +  Ii.5l99998x  +  .0Ui43000. 


Example  8 

f(x)  =  x^  -  .9Ux^  +  .6x^  +  2.99X  +  10.1i5 
(from  Hitchcock's  example  using  Lin's  I9UI  iteration  [7^^]) 


r.'s  are ; 

1 


(y.  '  s  are: 

1 


-l.Ol;  +  1.08i 
1.51  +  1.55i 

-  .076li8Iii; 
.0391521 

-  .OI19OOOO 
0.0 

0.0 


\,\ 


\\ 


-.00285U9  +  .00879511 
-.00281;88.+  .00833801 


,00285ii9  +  .OOi;98l5i 
,0028Ii89  +  .001499211 


Therefore  the  zeros  are:       , 

-1.0ii28ii88  +  1.07116621 
1.5l28i;89  +  1.51;50079i 

Check:     n(x-r?)    =  x^  -   .9U00002x^  +   .6000012x^  +  2.989999x 

+  10.U500015. 


hS 


Example  9 


f(x)   =  x^  -  5.92x^  +  13.66lx^  -  I5.i;05x  +  7.5308512 

estimated  r    's:  1+i 

27  2 

a   's  are:  -.ii69lii88 

.5950000 
-.3390000 
- . 0800000 
0. 
0. 


.0085372  T  .019251  .0021^258  .09707li5 

.0099651  ~  .01997171       .002l;9U6  .099777ii 

.0099967  +   .01999861       .0025000  .0999938 

.0100000  +  .02000001  .0025002    .0999996 

Therefore  the  zeros  are: 

1.0099999  +  .98000001  as  compared  to  the  1.01  +  .98! 
1.9500002  +  .OOOii553i  actual  zeros:      1.95;  1.95. 

Example  10 

.   f(x)  =  x^  +  256. 8x^  +  151. 9x^  +  256. 5x^  +  88. 5x^  +  h.^x  +  .05 
(from  Example  h   of  Chapter  II) 


r    ' s  are : 

-.U8U625  +  1.095U 

-.0157900" 

-.0U28000 

.iil486598 

-256.2111^02 

0!.  '  s  are: 

.I6IU605 

1 

13.990i;507 

238.3653011 

-5.0881a71 

2.1765129 

0.0096817 

0. 

0. 

1^6 


.2178881  + 

.  20095ii8i 

.OOOOOOii 

.0000583  - 

.lili55926 

. 0000761 

.3283h07  + 

.2092l5Ui 

.0000152 

.  00102l;5  - 

.7035223 

. 0000761 

.3809192  + 

.16606231 

.0000013 

- . 0002359  - 

.8032323 

. 0000761 

.3872277  + 

.  1381^07  81 

. 0000008 

-.00011;77  - 

.7891522 

. 0000761 

.3795185  ~ 

.Iiil3l487i 

-.0000008 

-.00011461  - 

. 7598811 

.0000761 

.3755183  " 

.lii970398 

-.0000009 

-.0001628  - 

.7581017 

. 0000761 

.3770083  + 

.I50ii719i 

-.0000009 

- . 0001667  - 

.7652827 

. 0000761 

.3783201  ~ 

.lU85263i 

-.0000009 

-.0001627   - 

.7672360 

. 0000761 

.3782191  + 

.lU92l4iiOi 

-.0000009 

-.0001609  - 

.7658265 

.0000761 

.3778706  + 

.114829301 

-.0000009 

- . 0001615  - 

.7650339 

. 0000761 

.3778220  + 

.Iii852li2i 

-.0000009 

-.0001621  - 

.7652305 

.0000761 

.3779017  ~ 

.11829301 

- . 0000009 

-.0001620  - 

.765i;73U 

. 0000761 

.3779301;  + 

.Iii81i091i 

- , 0000009 

- . 0001619  - 

.765U73U 

.0000761 

Therefore  the  zeros  r 

'.  are: 

1 

-.10669ii5 

+  .9ii69909i 

-.0157909 

-.Ol;296l9 

-.3168139 

- 

256. 21106  3li. 

Check:     i^Cx-r?)    =  x^  +  256.8000l83x^  +  I5l.90li57l5x     +  258. 50607  30x^ 
+  88.5l5967Ux^  +  I4.5009227X  +  .0500106. 
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12x^  +  50x^ 


Example  11 

f(x: 

(from  Example  5  of  Chapter  II 


r.  ' s   are: 

1 


en.  '  s  are: 

1 


70x- 


lilx     +  2x  +  390 


-1.10688142  +  1.2331 
2.977IO6I4  +  2.09291 
3.2595779  ~ 
14.9999778 

■  202.58031146 

90.920I4O35 

-15.9860919 

1.1802006 

O.OOOII483 

0.0 

0.0 

0.0 


hi 


\.\ 

V\ 

\ 

\ 

0959058  ~ 

.19596831 

.01iil5Ul;  + 

.08650li 

-.220ll;37 

.0000232 

105l9li8  T 

.227i;607i 

.0216165  + 

.09238U6i 

-.255I1O66 

.0000226 

1067177  + 

.2322ii30i 

.0227807   + 

.09281671 

-.2590590 

.  00002  2I1 

1068525  ~ 

.2  3289591 

.0228781  T 

.09289071 

-.2597132 

.0000221; 

1068811  ~ 

. 23298561 

.0228919  ~ 

.09289831 

-.2595686 

.  00002  2ii 

Therefore  the  zeros  r " ' s  are: 


-1.0000031  +  l.OOOOlliU 
2.9999982  +   2.00000171 
3.0000092 
5.0000002 


as  compared  to  the   actual  zeros: 

-1+1;  3+21; 


3; 


C.  As  a  Method,  Independent  of  Other  Methods, 
for  Finding  the  Zeros  of  a  Polynomial 


Example  12 

f(jc)    =  x^  -  5.0000002x^  +   3.0000029x^  -   2.000022ipc  +  1.0000006 

Noting  that  f(0)  >  0  and  f(l)   <  0  and  that  f(ii)   <  0  while 

f(5)  >  0,    so  that  there   are  real   zeros  between  0  and  1  and  between 

h  and  5,   we  take  for  Initial  guesses     .5  and  h.S  i'or  the  real  zeros  of 

f(x),    and  since  f (x)/(x-.5)(x-i4.5)  ~  x     +   .75  whose   zeros  are  ^  +   .871, 

we  take  +    .871  as  the  remaining  values  of  r. 's. 

Here  Q-    '  s  are:  -.70302iiii 

1.78Ui;775 
-.0068971 
- . 0000002 

.0 

.0 


Ii8 


\ 


Tl, 


v\ 


,0855367 
.0886ii85 
,08903Ut 


.019U928  +  .21180311 

.OIOI4368  +  .2i|06905i 

.0069015  ~  .2ii23776i 

.00683U9  ~  .2i|22557i 

-.0889739  .00687li7  ~  .2^225161 

Therefore,   the  zeros  are:  .57522ii.9 

i;.i;11026l 
.  00687  U7   +   .6277l|8i;i 


.0165513 
.0721217 
.0751^015 
.0752518 
.07522U9 


as  compared  with  the   actual  zeros 


.5752290 
U.UIO253 
.0068730  +   .6277^601. 


Example  13 

h(x)  =  x^  +  8.0000038x^  +  8.0799952X  -  15.8399957 
Considering  xV8. 0000038  +  (x^  +  1.0099990x  -  1.9799985), 

the  zeros  of  the  quadratic  part  are  readily  obtainable,  in  fact, 

2 
~  -2  and  1  since  the  quadratic  polynomial  ~  x  +  x  -  2.   So  tempor- 
arily truncating  the  cubic  polynomial,  we  get  for  our  input 


r^-. 

3:          -2  and  1 

a    '5 

3:           ,0200015 

X 

.0099990 

.0 

.12U9999 

\ 

\ 

^o 

1 

2 

3 

-.3333319 

-.0516668 

I.38I9987 

-.5879U87 

- .  Oii885ll 

1.587lii8l 

-.61U2955 

-.0I466220 

I.66129I17 

-.62U3220 

-.0I167U98 

1.6690061 

-.625ii606 

-.Oi;66695 

1.6721U29 

-.6258970 

-.Oii667W 

I.672I4822 

-.6259U72 

-.0U6671I4 

1.6726190 

-.6259662 

-.01^66716 

1.6726339 

h9 


Therefore  the  zeros  of  h(x)  are:    -2.6259662 

.953328U 
-6.3273725 

as  compared  with  the  actual  zeros:      -2.6259698 

0.9533286 
-6.3273627. 


Example  llj 


h(x)   =  x^  +  9.9999998x^  +  .00000l5x^  -  9. 999991 Sx^ 

+   .00000ii8x  -  19.9999990 
The  zeros  of  h(x)    are  the  same  as  those  of  g(x)   =   .Ix^  +  f(x) 


where  f(x)  =  x^  +  .0000002x^ 


,9999998x  +  .0000005x  -  1.999999 
^  -  x^  -  2 


1     0 
Since  the  zeros  of  x  -  x  -  2  are  readily  found,  we  use  the  technique 

of  temporarily  truncating  g(x)  to  obtain  the  zeros  +  v^and  +  i  of  f(x) 


So  take  r   ' s: 

+ 

i 

+ 

l.iilii21 

a.  's: 

1 

-.0000100 

. 0000005 

-.0000099 
.0000002 
.0 

.1 

11    ,   11 
1'     '2 

\ 

\ 

^5 

.0166668 
.0167226  + 
.0165776  ~ 
.0165772  ~ 

.001080i;i 

. 00105551 
. 00105611 

-0.0666629 

-  .0581;762 
-..0582909 

-  .0582852 

-.0666699 
-.0792109 
-.0790671 
-.0790720 

.0999992 
.IOU12U2 
.10i|2036 
.IOU2031 

Therefore,   the 

zeros  are 

.0165772  + 
1. 35592^7 
-l.ii932820 
-9.8957968 

.9989ii39i 

as  compared  with  the   actual 

.  zeros: 

.0165771  + 
1.35592ii9 
-I.li9328l9 
-9.8957971. 

.9989li39i 
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Conclusions 

¥e  have  thus  established  that  the  iterative  formulas  proposed 

in  this  paper  converge  v^enever  the  fee  |' s  in 

^n   n+1        ,  , 
n+1    V-  /     v  1-1 

n+2        1  ^  i  i^ 

are  sufficiently  small.  As  a  means  of  solving  for  the  zeros  of 

n+1 

E  (a.+«. )x^~^  , 
111^     ' 

that  isj  Q-  ,p  =  Oj  that  condition  amounts  to  requiring  that  the  esti- 
mates  r. 's  of  the  actual  zeros  r. 's,  vrfiich  are  input  of  the  program,  be 
not  too  far  off.  However,  as  shown  in  Examples  10  and  11,  the  process 
frequently  converges  even  vri:en  the  a   's  are  as  big  as  a  three-digit 
figure. 

The  iteration  is  then  particularly  effective  as  a  corrective 
measure  to  other  methods  which  obtain  the  zeros  one  at  a  time  and  thus 
are  "contaminating"  procedures.  Note  from  the  example  that  as  such, 
this  noncontaminating  iterative  process  usually  obtains  the  actual 
zeros  after  only  a  few  iterations. 

As  a  method,  by  itself,  of  obtaining  the  zeros  of  a  polynomial, 
it  is  especially  applicable  to  a  polynomial  of  degree  3  or  one  of  higher 
degree,  wherein  ignoring  the  leading  term  gives  a  polynomial  of  lower 
degree  vrfiose  zeros  are  readily  obtainable  or  are  known.   (Examples  13 
and  lii  illustrate  this  point.)  Then,  after  the  zeros  of  the  truncated 
polynomial  are  found,  the  zeros  of  the  original  polynomial  can  be 
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determined.  Here,  (x     .  /  0     and  the  iterative  process  seems  to 
converge  when  Jcc  ^pl  ^^  quite  small,  that  is,  when  the  second  coeffi- 
cient of  the  original  monic  polynomial  is  rather  large  in  magnitude. 
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APPENDIX  A.  MACHINE  PRGGRAM3  IN  FORTRAN  II  FCR  IBM  709 
LIN'S  PENULTIMATE  REMAINDERING  WITH  TRANSLAFION  FOR  QU 

ADRATIC  OR  L 
IINEAR  FACTOR, WIIH  GCD  ALGC  TO  DETERMINE  PRESENCE  OF  MU 

LTIPLL  ZEKCS 
Al{20)  ,  6(2,20)  ,  C(20) 
Y(2) 
10,  N ,  J  L  I  M 


14 


A9 


31 


79 


53 


57 


55 


61 
59 


63 


67 


65 


0(20)  ,  AD(^iC) 


14,  (A(  I  )  ,  1  =  1, N) ,  H 


I0,\  OF  A'S  BY  AD'S 


LIMENSION  A(20), 

DIMLNSIGN  G( 10) , 

REAUINPGTTAPES, 

FORMAT (2 12) 

REACINPLTTAPE5, 

F0RMAT(4F15.7) 

GCD  ALGORITHM  TO  UEFERMINE  HREShNCE  CF  MULTIPLE  ZERCS 

CO  49  1=1, N 

Al(  I  )=  A(  I) 

Nl  =  N-1 

AN=  N 

DO  51  1=1, Nl 

AI  =  I 

AD(  I  )=  AI*A1(  1  +  i )/AN 

NN  =  N 

M  =  N1 

PERFORM  DIVI 

Nl=  NN-1 

DO  53  1=2, M 

MN=  NN-M+I 

C(MN)=  0.0 

MN=NN-M+1 

C(MN)=  1.0 

CO  55  I=M,N1 

DELTA  =  O.C 

CO  57  K=1,M 

MN=  NN+M+1-I-K 

CELTA=  C(MN)*AD(K) 

NI=  NN-I 

NIM=  NI+M 

C(NI  )=  AKNIM)-  DELTA 

NM1=  NN-M+1 

C(NM1)=  1.0 

DO  59  1=1, M 

DELTA  =0.0 

DO  61  K=1,I 

IK=  I+l-K 

DELTA=  CELTA+ 

D(  I  )  =  Aid  )- 

CHECK  WHETHER 


+  DELTA 


AD(K)»C{ IK) 

DELTA 

DIVISION  NEED  BE 


REPEATED 
INED 


OR  bCD  IS  CBTA 


DO  63  1=1, M 

MM=  M-I+1 

IF( ABSF(C(MM) )-1.0E-5)  63,63,65 

CONTINUE 

WRITFCUTPUTTAPE6,  67,  (ACd),  1  =  1, M) 

FORMAT  (14H  THE  GCD  IS    /{6Fi2.6)) 

GO  TO  81 

IF(MM-l)  69,69,71 


^h 


69    l^RITfE0UTPUTTAPE6,     73 
73    FCRN!AT(gH    GCC    IS    i) 

GO    TC    91 
71    l\iN  =  K 

M=    MM-i 
DO    75    1=1, NN 
75    AU  I)-    AD(  I  ) 

CO    77    1=1, M 
7  7    AD(  I  )=    C(I)/D(f^M) 

GO    TC    79 
81    DO    83    1=2, M 

MN=    N-l^  +  l 
83    C(MN)=    0.0 
N"N=    N-M+1 
C(Ml\)=    1.0 
Nl=    N-1 
CO    85    I=M,M1 
CELTA=    0.0 
DO    97    K=1,M 
MN=    N+M+1-I-K 
97    DELTA=    C(MN)«AD(K)     +    DELTA 
NI=    N-l 
NIM=NI+M 
05    C{NI)=    A(NI^')-DbLTA 
N  =  H-H 

CO    89    1=1, N 
89    A{I)=    C(I) 

91  WRITE0UrPUTTAPE6,  32,  (A(I),  1=1, N) 
32  F0RKAr(46H  THE  CGEFF  OF  POLYhCMIAL 

IF(N-2)  1,1,93 
PERFCR^^  TRANSLAIION  ON  X  OF  F{X) 
93  CO  174  1=1, K 

174  C(I  )  =  A( I) 
DO  175  1=  1,N 
ND  =  N+2-I 
D{NC)  =  1.0 
CO  176  K=I,N 
NK=  N+l-K 

176  D{i\K)  =  C(NK)  +  H*0(NK+iJ 
A 1  n  )  =  D  (  1  ) 
IF(I-r^)  2CC,  175,175 

200  NI=    N~I 

CO    2Ci    K  =  1,M 

201  C{K)    =    D(K+1) 

175  CCNTI;<iUE 
A1(N+1)=    1.0 

WRITE0UTPUTTAPE6,  32,  (Aid),  I-1,N) 
FOR^'  THE  1ST  TRIAL  DIVISOR 
IF{XM0DF(N,2))  22,23,22 

22  M  =  1 

GO  TG  25 

23  F  =  2 

25  IF(A8SF(Al(F+l))-i.0E-5)  27,27,29 


3EING  FACTORED 
/(4F15.7)  ) 


ARE 


55 


27 

20 

29 
31 


34 
16 


36 


40 

38 

44 

50 
46 

48 
42 
56 

64 
62 

54 
66 

68 

106 

70 

74 

82 


I   87 


WRITtOU 
FORMAK 
GO  TO  1 

00  31  I 
B{  I,l)  = 
PERFCRM 
J=l 

WRITfcOU 
FORNAK 
WRITEOU 
DO  36  I 
MN  =  N- 
C(MN)= 

C(MN)= 
N1=N-1 
DO  38  I 
DELTA  = 
DO  40  K 
MN=  NtM 
DELTA  = 
M  =  H- 
NIM  ^  N 
C(NI  )  =  A 
IF( J-1) 
T=AeSF( 
IF( T-1. 
IF(M-l) 
T=  ABSF 
IF( T-1. 
IF( J-JL 
IF( AESF 
B( 1, J+1 

1  F  { ^'-  1  ) 
B(2,J+1 
J  =  J  +  1 
GO  TC  3 
WRITEOU 
FORMAT( 
GO  TC  1 
OBTAIN 
X=  -B( 1 
WRI  TEOU 
FORf^AT 
GO  TC  7 
CO  74  I 
G(  I  )  = 
G(  I)  = 
CALL  QU 
hH  =  (0 
hH  =  H 
DC  87  I 
Y(  I  )  = 
WRITEOU 


TPUTTAPE6,  20 

3gH  G  DIVISOR,  CHANGE  STARTING  VALUES  OF  B) 

=  1,M 
Al(  I  )/Al(K+i) 

DIVISIOiM 

TPUTTAPEo,  16,  J 

IX,  18H  ITERATION  NG.  IS  ,12) 

TPUTTAPE6,  14,  (C(I,J),  1  =  1, M) 

=  2,M 

M  +  I 

O.C 

+  1 

1.0 

=  N,N1 

0.0 
=  1,M 
+1-I-K 

C(MN)»B(K, J)+DELTA 
I 

-I+M 
1(NIM)-DlLTA 

42,42,44 
B(l,  J)-B{  1,J-1) ) 
OE-3)  50,50,46 

1,68,46 
(B(2,  J)-u(2,J-i)  ) 
OE-3)  7C,7C,48 
IM)  42,42,54 
(C(  1)  )-1.0E-5)  27,27,56 
)  =  AK  1)/C{  I) 

1,  62,  64 
)=  {A1(2)-C(2)*B( I, J) )/L(l  ) 


rPUTTAPE6,  66, N 

40H  ITERATION  NOT  CONVERGING  FOR  DEGREE     /I2J 

THE  ZEROS 

,J)+  H 

TPUTTAPE6,  106,  X 

(  13H  ZEROS  ARE    /(2F12.5) ) 

2 

=  1,2 

(0.0,0.0) 

B{  I,J) 

ADRA(G,Y) 

.0,0.0) 

=  1,2 

Y I  I )  +  HH 

TPUTTAPE6,  106,  Y(l),  Y(3),  Y(2),  Y{4) 


% 


72  H=H-^ 

IF(N)  Itl.a-!. 

84  CO  80  1  =  1, l>J 
80  Aid)  =  C(  1  ) 

IF(N-2)  I,  150,  12 

150  H=2 

DO  151  1=1,2 
G{I)=  (0.0,0.0) 

151  G (  1  )  =  A  1  (  I  ) 
GO  TC  82 

12  WRITtOUTPUTTAPiE6,  32,  (AUI),  1  =  1, N) 
GO  TO  2  3 
END 

SUBRCUriNE  gUACRAlCY) 
LI  WENS  ION  G{  10),  Y(2) 

DELTA=  SQRTF(G(2)«»2-(4.0,0.0)*G( 1) ) 
Y(l)  =  (-G{2)+D£LTA)/{2. 0,0.0) 
Y(2)=  (-G{2)-DELTA)/(2. 0,0.0) 
RETURW 
END 


NDFILE 


For  the  data,  n  =  the  degree  of  the  polynomial  f(x)  (<  20) 

jlim  =  upperbound  for  the  iteration  number 

h  =  amount  of  translation  on  x 

n+1 
{A(i)]  =  real  coefficients  of  f(x)  =  S  A(i)x^~ 

1 

If  GOD  =  1,  Penultimate  Remaindering  is  applied  to  f (x+h) . 

If  GOD  is  a  polynomial  of  positive  degree,  Penultimate  Remaindering 
is  applied  to  f/GCD,  in  order  to  have  simple  zeros.  The 
multiplicity  of  the  zeros  found  can  later  be  determined 
either  by  studying  the  GOD  or  by  repeatedly  depressing  f 
by  the  zeros  obtained. 


57 


lA 
CI 
I      CI 


L,G( iC), 


A  NCNCONTAMINATING  ITERATIVE  PROCESS  FOR  DETERN'.IMNG  T 

HE  ZEROS  OF 
POLYNOMIAL 
PENSION  NT{5) 

MENSION  ACIO),  ALPHAdZ),  C(IO),  00(10,10,20),  D(10) 

,  ETA(11,^0) 
H{10,  5),  R(10),  RR(ll),  X(ll),  Y{2),  PNUM(IG) 

,  PNOFKlL,  5 
(10),  SUM(5,10),  SSLN(5) 
TTAPE5,  1,  N,  NS,L 
313) 


TTAPE5,  2,  (A(I),  A(I+10),  1=1, N) 

TTAPE5,  2,  (R(I),  R(I+1U),  1=1, N) 

TTAPEti,  2,  ALPHA(N2),  ALPHA{N2+12) 

POTTAPEO,  7 

6H  HERE  ARE  THE  A'S,  R'S  +  ALPHA'S,  RESPECTIVE 

LY) 
PUTTAPE6,  2,  (A(I),  Ad  +  lO),  1  =  1, i\) 
PUTTAPEo,  2,  (R(I),  R(I+10),  1=1, N) 
F15.7) 
3,4 

TrAPE5,  1,  (NT(I),  1=1, L) 
LPHA'S 
R(  1) 
2,N 
(  I) 

TIP(0,D,G,I-1,1 J 
1,1 
J) 
1,N 

=  A{J)-0(J) 
)=  (CO, 0.0) 

PUTTAPE6,  2,  (ALPHA(J),  ALPHA(J+12),  J=1,N2) 
LE  R«S  AND  R(N+1) 

)  =  GIN)  -  ALPHA(N2)»(G(N-1)-G(.\)*»2)  +  ALPHA! 

N) 
7,27,12 
1,NS 

=  (0.0,0.0) 
=  (1.0,0.0) 

1,N2 
=  PNUM(K)-  ALPHA( I )»R(K)*»( I-l) 
1,N 
16,15,16 

=  CENOM(K)»(R(K)-R{ I ) ) 

1)=  PNUM(K) /DENOM(K) 

ALPHA(N2) )  )  19,  19,20 

(0.0,0.0 

(1.0,0.0) 


2),  CENOM 

100 

REACINPO 

1 

FORMAT  ( 
Nl  =  N+1 
N2  =  N+2 
REACINPU 
REACINPU 
REACINPU 
URITEGUr 

7 

F0RNAT(4 

V^RITEOUT 
WRITECUT 

2 

FGRMAT(4 
IF(L)  3, 

4 

REACINPU 

C 

OBTAIN  A 

I 

3 

C(  1)  =  -1 

CG  10  1= 

I 

C( 1)=  -R 

I 

CALL  ^'UL 
CC  10  J= 

I 

10 

C{J)=  G( 
CO  11  J= 

I 

11 

ALPHA( J) 

I 

ALPHA(N1 
WRITtCUIi 

C 

FOR  SIMPI 
J=l 

I 

ETA(M,2 
IF(NS)  2 

12 

CO  13  K= 

I 

PKUM(K) 

I 

DENCM(K) 
CO  14  1= 

I 

14 

PNUI^IK) 
00  15  1= 
IF  (K-I) 

I 

16 

CENC(^(K) 

15 

CONTINUE 

I 

13 

ETA(K, J+ 
GO  TO  17 

I 

18 

IF( AeSF( 

I 

20 

PNUM2  = 

I 

PNUN3  = 

58 


21 


23 
24 

22 
26 
25 
19 
28 


I   31 

30 

33 
I  34 
I  32 
I   35 


I 


29 
17 
27 
37 

36 
38 


39 

41 


I 

44 
I  220 


47 

46 

48 
45 
43 

49 


CO  21  I 
IN  =  H- 
PNUf^2  = 
Pt\U^3  = 
PNU!^4  = 
CENCM 
IP(NS) 
CG  24  I 
CENGMl 
IF  (L) 
DO  26  I 
DEKCKl 
ETACNl, 
IF(NS) 
CO  29  K 
CEMCN(K 
CO  30  I 
IF(K-I  ) 
CENOFCK 
CQNl  INO 
IF(L)  3 
CO  34  I 
CENGMCK 
IF( ABSF 
DENDM(K 

ETA(K, J 
WRI lEOU 
IF(ABSF 
WRITECU 
WRIT  ECU 
FGRI^AT 
FOR  I^UL 
IF(L)  3 
IF( J-1) 
NSS  =  N 
NSl  =  N 
CO  4  2  I 
MI  =  N 
NSI  =  N 
C(  1)  = 
IF(NS) 
CO  45  K 
IF( J-1) 
D( 1)=  - 

GO  rc  4 

C( 1)=  - 
CALL  MU 
DO  48  K 
C(KK)= 

CONTINU 
NCEG=  N 
IF { J-1) 
CO  51  K 


=  1,N 
I  +  l 

PNUK2  +  ALPHA{  I,\)*(-ALPHA{N2)  )*»l 

((1.0, CO)  +  R(  n»ALPHA(N2)  )*P!\ILf3 

(i.CO.G)  -  PNUM2-PNUM3 
=  ALPHA(N2) 
22,22,23 
=  1,NS 

=  DENOMl  *(  ( i.O,  .G)+ALPHA(N2)«(R(I )+cTA( I ,J)  ) ) 
25,25,22 
=  1,L 

=  U£K0M1*SSUM(I) 
J+l)=  PNUM4/0EN0M1 
27,27,28 

=  1,NS 
)  =  (  1.0,0.0) 
=  1 ,  IN  S 

31,30,31 
)=  DENQM(K)»(R(K)-R( i)-ETA(I,J)) 
£ 

2,32,33 
=  1,L 

)  =  CEN0K(K)»SUM(1,K) 
(ALPhA{N2)))  29,29,35 
)  =  DcNOM(K)  »(  (  i.O,0.0)+ALPHA(!\i^)*(K(K)-tTA(M 

,  J)  )  } 
+  1)  =  PNGM(K)/DEiNiUK(K) 
TPUT1APE6,  2,  (ETA(K,J+1), 
(ALPHA(N^)))  36,36,3? 
TPUTTAPE6,  2,  ETA(N1,J+1), 
TPUTTAPE6,  38 

(32H  ABOVE  ARE  EIAS  FOR  SIMPLE 
TIPLE  ZEROS 
9,39,41 

94,94,40 
S  +  1 
S  +  1 
=  l,L 
T(I) 
SS 

(  1.0,0.0) 
43,43,44 
=  1,NS 

220,220,47 
R(K) 
6 

R(K)  -  EIA(K,J) 
LTIP{C,C,G,K-1,1) 
K  =  1,K 
G(KK) 
E 
S 

49,49,50 
=NS1,N 


ETA(K,J+2i),  K=1,NS) 

ETA(N1,  J  +  21) 
ZlRLS) 


59 


53 
I   52 

I 


I   54 
51 


1   61 


50 


57 

I   58 
I 


59 
56 


67 


I  210 
I  55 
I 


9 
63 


I 
I 
I   64 


65 


66 


IF(R(K)- 
IF(R(K+1 
C( 1)  =  - 
CALL  MUL 
NCEG=  ni 
DO  54  KK 
C(KK)=  G 
CONTINUE 
LO  61  K= 
PNUMUK, 
Kl  =  K-1 
DO  61  M= 
Ml  =  M-1 
PNUN*!  (K, 

GO  TC  55 
CO  56  II 
NTH  =  N 
IF(  I-II  ) 
CO  58  K= 
C(K)=  h( 
CALL  MUL 
NiDEG  =  N 
CO  59  K 
C(K)  =  G 
CONTINUE 
IF( ABSF( 
CCN=  (1. 
GO  TO  55 
CON  =  (  1 
RT  =  R(N 
CALL  DER 
CO  60  K= 
Kl  =  K-1 
CENCM2  = 
PNUM2  = 
PNUM3  = 
IF(Kl)  9 
IF( J-1) 
CO  64  M= 
K1M=  Ki- 
KM=  K-M 
RT=  R(NS 
CALL  DER 
PNUM3  = 
IF( J-1) 
CO  66  M= 
KIM  =  Kl 
KM=  K-M 
PNUN2  = 
AK=  Kl 
PNU^Z  = 
PNUI^'B  = 
CEN0M2  = 


R(NSI))  52,53,52 

0)-R(NSi+10) )  52,51,52 

R(K) 

TIP(C,D,G,NDtG,l) 

EG+1 

=1,NCEG 

(KK) 

1,NTI 
I)  =  (0.0,0.0) 

K,N2 

I)  =  PNUMKK,  I  )+PROD{Ml,Kl)»ALPHA(M)»R(NSI  )*»' 

M-K) 

=  1,L 
T(  II  ) 

57,56,57 
1,NTII 
K,  II  ) 

TIP(C,D,G,NDE&,NTII ) 
CEG  +  Nl II 
=  1,NDEG 
(K) 

ALPHA(N^)))  210,210,67 

0,0.0)  +  ALPHA(N2)»(R(NSI  )-ETA(.M,  J)  ) 

.0,0.0) 
SI) 

IV(C,C,X,N-NTI,  RT) 
1,NTI 

X(  1)»PK0C(K1,K1) 
(O.CO.C) 
(O.CO.O) 
,9,63 
60,60,8 
1,K1 
M 

I  ) 

IV(M,C,X,N-NTI,  RT) 

PNUM3-PR0D(Kl,KlM)»X(Ni+l  )»C0(KM,I,J+1) 

60,60,65 

1,K1 

-K 

PNUM2  +  PR0U(K-2,  K  itM  )  »C0  {  KM  ,  I  ,  J+ 1  )  «-X  (  H  ) 

-ALPHA  (rM2)»AK*PNLM2 
CCN»PNUM3 
CON  »DtNaM2 


60 


I       60    CC(K,I,J+1)     =     (FNUM1(K,I)     +PIMLM2  +  PNU^3 ) /0£:NGNi2 

WRITEOUTPUTlAPEo,     2,     ( CG ( K , I  ,  J+ i )  ,CG ( K  ,  i  ,  J  +  2 1 )  ,    K=l,Nr 

1  ) 
WRITEGUrPUTTAPE6,     6'9,     I 


271 


I 

I 

I 
I 
I 

I   73 


69  FCRf'ATOSH  ABOVE  ARE 
CATA  NEEDED  FOR  ,'NltXT 

lOS 
IF(^SJ  62,62,70 
DO  71  K=l,NS 
SUH(I,K)=  (0.0,0.0) 
DO  271  KL=1,NTI 
SUM(I,K)=  SUM(  I  ,K)  + 


C'S  FOR  f^LLTIPLE  ZERO    /I3) 
ITERATION  OF  BOTH  THE  SIMPLE  AND 
MULTIPLE  ZER 


70 


C0(KL,I,J+1)»(R(K)-R(NSI))»«(KL-1 

) 
(R(K)-R(NSI  )  )»»KTI 


71  SUI^(  I,K)=  SUM(  I,K) 
62  IF( A3SF{ALPFA(N2)  )  )6H,6&,72 

72  SSOf(I)=  (0.0,0.0) 
CO  73  K=i,NTI 
SSUM(I)=  SSUMd)  +  C0(K,I,J  +  1)»(  (-ALPHA(N2)  )»*(MI-K  +  i 

)  )  «  (  (  i  .  0  ,  0  .  0 
1)  +  ALPhA(i\2  )*R(NSI)  )*»(K-1  ) 

SSUH{I)=  SSUMd)  +  ((1.0,0.0)  +  ALPHA(N2)»R(NSI)  )»»i\TI 

CO  74  K=1,NTI 

C(K)  =  CG{K, I,J+1) 

LG  75  M=  1,NTI 

ND  =  'MTI+2-M 

C(NL)  =  (1.0,0.0) 

DO  76  K=  iy,NTI 

NTK=  NTI+l-K 

D(NTK)=  C(NTK)-R(NSI )«D(KTK+1) 

H  (  M  ,  I  )  -  C  (  1  ) 

IF  (N-NTI)  200,75,75 

NTH  =  iMTI-M 

DO  201  K=i,NTM 

C(K)  -  C(K+1) 

CONTINUE 

MSS=  NSS+NTI 

TEST  FOR  CCNVtRGENCc 

IF(J-l)  94,gA,7d 
78  DO  80  1=1, L 

NTI  =  NT( I ) 

DO  80  K^l.NTl 

T  =  A8SF(C0(K,  I,J  +  1)-C0(K,I  ,  J)) 

IF(T-l.GE-4)  80,80,82 
30  CONTINUE 

IF{NS)  81,81,40 

CO  83  K=1,NS 

T  =  ABSF(ETA(K,J  +  1)-ETA(K,J)  ) 

IF(T-1.0E-4)  83,83,82 
CONTINUE 

IF( A8SF(ALPhA(N2) )  )  84,84,77 

T  =  ABSF(ETA(IM1,  J+1)-ETA(N1,  J)) 

IF(T-1.0E-4)  34,84,82 
84  WRITL0UTPUTTAPE6,  85 


63 
74 


I   76 


200 

201 
75 
42 


40 


83 
61 
77 


61 


85 

87 
I 

88 
I  86 
I   99 

97 

89 


9L 
I  92 
I 


I 


93 

90 

82 
94 
96 
79 

95 
98 


I  101 

I  102 

I 

I 

I 


I  103 
I 


106 
107 


FORMAT 
IF(NS) 
DO  88 
RR(K)= 
WRITEO 
IFCAOS 
RR(Nl) 
WRITEO 
IF(L) 
NSS=  N 
CO  90 
NSI  = 
IF(NT( 
CO  92 
G(K)  = 
CALL  Q 
CO  93 
RR{K)= 
WRITEO 
NSS  = 
GO  TC 
IF( J-1 
WRITEO 
FOR^'AT 
J  =  J  + 
GO  TG 
WRITEO 
FORMAT 
GO  TC 
ENC 

SUBRGU 
CIMENS 
NP2  = 
NP2  = 
NM=  NP 
CO  ICl 
C(K)  = 
CO  1C2 
D(K>  = 
C(NP+1 
C(MP+1 
CC  1C3 
G{  I  )  = 
CO  ICi 

IJ=  I- 

G(I  )  = 

RETURN 
END 

FUNCTI 
MPRCC 
IF(NU) 
IFCNV) 
DO  108 
NUI  = 


UH 

6,86 

1,NS 

{K)  + 

POTT 

ALPh 

-(1 
PUTT 
0,10 
1 

1,L 
S 

-2) 
1,2 
C(K, 
CRA( 

1,2 
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108  KPRCC  =  NPRCD»NLI 
I  105  PRCC  =  NPROD 
I      RETURN 

END 

SUBROUTINE  DER  I V ( M, G , X , ME » R T ) 
I      DIN!ENSION  X(  11),  G(  10) 
I      X(M+1)  =(O.CfO.O) 

IF(i^-ME)  110,110,111 
110  Ml  =  M+1 

DC  112  I=M1,  ME 
i  il2  X(M1)  =  X(M1)+  PRGDC I-1,M)«G( I)*RT««(  I-Ml) 
I      X(M1)  =  X(M1)  +  PRaD(ME,  M ) *RT»« ( ME-M ) 
I  111  RETURN 

END 

SUBROUTINE  GUACKA(G,Y) 
I      DIMENSION  G{  10) ,  Y(2) 

I      DELTA  =  SeRTF{G(2)»*2-('^.0,G.C)*G{  1)  ) 
I      Y(l)  =  (-G(2  )+DtLTA)/(2. 0,0.0) 
I      Y(2)  =  (-G(2)-CLLTA)/(2.C,C.C) 
I      RETURN 

END 


NDFILE 


For  the  data,  n  =  degree  of  the  polynomial  f (x)  (<  lO) 

ns  =  number  of  simple  zeros  of  t(x) 

L  =  number  of  multiple  zeros  of  t(x) 

n+1      . 
{A(i)}  =  coefficients  of  f(x)  =  E  A(i)x^" 

{r(i)}  =  zeros  of  t(x)  with  the  simple  zeros  listed  first 
and  the  multiple  zeros  repeated  according  to 
respective  multiplicities 

[nt(i)}=  multiplicities  of  the  multiple  zeros  of  t(x) 
alpha  (n2)   =  coefficient  of  x^"^  term  of  g(x)  , 


APPENDIX  B 
COMMON  METHODS  FOR  DETERMINING  ZEROS  OF  A  POLYNOMIAL 

I.  For  Real  Zeros 

1.  Half  Interval  Search  [6] 

(a)  First  approximate  by  graphical  or  any  other  means  an 
interval  containing  a  real  zero  of  f(x),  say  (d,b). 
Assume  f(d)  <  0,  for  the  purpose  of  discussion. 

(b)  Test  vrtiether  f(— o-)  is  less  than  or  greater  than  or 
equal  to  0.  If  f(— 5-)  =  0,  then  the  required  zero  is 
obtained;  if  it  is  less  than  0,  then  the  required  zero 
is  in  the  interval  (— r-^  b),  whereas  if  greater  than  0, 
it  is  in  the  interval  {d,—r-) . 

(c)  So  if  still  necessary,  repeat  the  process  with  the  inter- 
val (d,b)  replaced  by  either  i-rr-,   b)  or  (d,  -p-)?  as 
the  case  may  be. 

(d)  By  continuing  the  process,  we  thus  narrow  down  the  inter- 
val as  much  as  we  want  to,  and  hence  locate  a  zero  to  an 
accuracy  of  any  reasonable  number  of  significant  digits. 

2.  Method  of  False  Position  [6] 

(a)  As  in  the  previous  method,  start  with  an  interval  (d,b) 
in  which  f(x)  changes  sign. 
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6U 


(b)  Linear  interpolation  is  then  used  to  find  an  approximate 
zero  X,  as  follows: 

Suppose  f(d)  <  0  and  f(b)  >  0.  The  line  joining  f(d) 
add  f(b)  intersects  the  real  axis  in  x^  where  x  is  deter- 
mined by  the  relation  f/KN  _  f(^)   "  t^^  ^^^^  ^^'   ^^ 

X  =  d  +  „r  ,\    j.7^\    •      (The  same  formula  is  derived  when 
f (dj-i(bj 

f(d)  >  0  and  f(b)  <  0.) 

(c)  To  determine  which  of  the  two  intervals  (d,x)  or  (x,b) 
contains  the  true  zero  x  ,  check  the  sign  of  f (x) : 

if  f(x)  >0,  then  x  is  in  (d,x);  if  <  0,  then  x^  is  in 
(x,b) J  if  =  0,  then  x  =  x. 

(d)  If  the  zero  is  not  yet  determined,  repeat  the  procedure 
with  (d,b)  replaced  by  either  (d,x)  or  (x,b) ,  as  the 
case  may  be. 

(e)  Continue  the  process  until  either  |f(x)|  <  eps  or  the 
length  of  the  interval  <  eps,  where  eps  is  a  predetermined 
small  positive  number. 

3.  Newton-Raphson  Method  [6] 

(a)  Start  with  an  arbitrary  approximate  to  the  zero,  say  x^  . 

(b)  To  improve  the  approximation,  draw  the  tangent  to  the 
curve  f(x),  and  determine  its  intersection  with  the  real 
axis,  say  x  .  The  relation  therefore  is 

f(x^)  f(x^) 

tan  a   =  f'(x^)  =  3^,  that  is,  x^  =  x^  -  j^^^y 
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(c)  The  iterative  sequence  can  now  be  set  up  as  follows: 

^n+1  \  ~   f  (x  ) 
n 

provided  f '  (x  )  -^  0.  . 
n 

(d)  Repeat  the  iteration  until  |f(x  -> )  I  *~  eps,  or  vdien 

|x    -  X  I  <  eps,  or  else  when  it  shows  signs  of  diver- 
gence such  as  when  |f'(x  -i  )  I  '^  eps  or  the  number  of 
iterations  exceeds  Jlim.  Here  Jlim  and  eps  are  predeter- 
mined values. 

II .  For  Real  and  Imaginary  Zeros 
1.   Penultimate  Remaindering  [1,2] 

(a)  Start  with  a  first  trial  divisor  to  f(x),  usually  the 
last  m  terms  of  fix),  reduced  to  a  monic  polynomial 
divisor,  vrfiere  m  is  the  degree  of  the  divisor  being 
sought . 

(b)  Divide  f(x)  by  the  approximate  divisor,  until  a  penulti- 
mate remainder  is  obtained,  that  is,  a  remainder  of 
degree  m. 

(c)  Taking  this  penultimate  remainder  divided  by  its  leading 
coefficient  as  the  next  trial  divisor,  the  procedure  is 
repeated. 

(d)  The  process  is  continued  until  the  number  of  iterations 
exceeds  a  certain  number  without  showing  sign  of 
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convergence,  or  else  when  the  coefficients  of  the  succeed- 
ing trial  divisors  converge. 

A  slight  variation  to  this  iteration  is  to  take  alternatively 
for  trial  divisor  the  penultimate  remainder  of  degree  m+1  and  the 
ultimate  remainder  which  is  degree  m,  each  time  reduced  to  monic  poly- 
nomial. 

2.  Lin' s  Ultimate  Remaindering  Iteration  [7,8] 

n+1    .  ^  m+1    .  ^ 

When  f(x)  =  2  a.x    is  divided  by  d(x)  =  E  b.x   ,  the 
11  1   1 

ultimate  remainder  has  coefficients  which  are  expressions  involving 

the  coefficients  a's,  b's,  and  the  quotient  coefficients  c's. 

If  d(x)  were  a  factor  of  f(x),  each  of  those  remainder  terms  is 

identically  0,  So  to  obtain  an  iterative  set  of  forraulas  for  the  b's 

lAich  hopefully  will  converge  to  the  exact  coefficients  of  d(x) ,  we 

set  each  of  the  terms  in  the  expression  for  the  remainder  to  0,  and 

solve  for  the  b's  as  functions  of  the  a's  and  the  c's. 

To  illustrate  for  the  case  when  ra=2,  the  quotient  when  f(x) 

is  divided  by  the  quadratic  d(x)  has  coefficients  expressed  as  a 

recursive  relation  c  ,  ^  =  a  ,  ^  -  b^c  ,  -  b^ c  ,  ,, 
n-k-1    n-k+1    2  n-k    1  n-k+1^ 

k  =  1,2,...  n-2  where  c  =0  and  c  ^  =  1:  v^ile  the  remainder 
'    '  n  n-1 

coefficients  are  r  =  a  -b„c^-b  c 
s  =  a-,  -b,  c^ . 

Solving  for  b,  and  h^,    after  setting  r  and  s  to  0,  we  get 

h  ^2"^l''2 

b^  =  —   and   b_  =  , 

1   c^         2     c^    ^ 
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from  which  the  iteration  formulas  are  derived,  namely, 

a2-b(l,j+l)c(2,J) 


^^^'^'^^  -  ^(1737 


b(2,j+l) 


c(l,J) 


For  the  first  set  of  values,  we  may  let  b(l,l)  and  b(2,l)  be 

—  and  — ,  respectively,  unless  a.^  =  0. 

^3      ^3 


3.  Luke  and  Ufford's  Extension  of  (1)  [9] 

Instead  of  deriving  the  formulas  that  will  enable  us  to  get 
an  estimate  of  one  factor  at  each  stage  of  iteration,  this  extension 
derives  the  formulas  that  will  display  the  estimates  of  two  or  more, 
possibly  all  the  factors,  at  each  stage  of  the  iteration.  This, 
therefore,  eliminates  the  contamination  of  one  factor  by  another  due 
to  division  in  the  depression  of  the  polynomials  and  perhaps  due  to 
the  inaccuracy  of  a  factor  obtained.  In  the  said  iteration,  if  some 
sets  of  coefficients  turn  out  to  be  convergent,  while  others  don't, 
then  depression  of  polynomials  will  be  inevitable,  so  the  rest  of  the 
factors  can  be  obtained  by  similar  procedure  from  the  depressed  poly- 
nomial. 

However,  to  attempt  to  write  a  general  expression  or  set  of 
formulas  to  separate  a  polynomial  into  any  number  of  factors  presents 
considerable  algebraic  difficulties.   So  the  iterative  formulas  are 
set  up  according  to  the  individual  problem. 


7    ^_^ 
To  illustrate  for  a  sectic  ^  a.x  ~  seeking  three  quadratic 

2  2 

factors,  let  the  quadratics  be  x  +  b  x  +  b  ,  x  +  d  x  +  d^,  and 

2  2 

X  +  e^x  +  e  .  Dividing  f(x)  by  x  +  bpX  +  b^ ,  we  obtain  as  in  Lin's 

Method  a  set  of  fonnulas  for  the  quotient  coefficients.  If  the 

^    i-1 
quotient  is  T,  c.x         with  remainder  r  x  +  r, ,  the  recursive  relations 


are: 


%^   %-   ^2 


C3  =  a^  -  c^b^  -  b^ 


c^  =  a^  -  c^b^  -  C3b^ 


^2  =  ^2  -  ^1^2  -  ^2^ 

^1  =  ^  -  '^l^ 

^    i-1      2 
Dividing  farther  E  c  x     by  x  +  d_x  +  du. ,   -we   derive 
1  ^  ^     i 

similar  set  of ^formulas  for  the  coefficients  of  the  quotients 

2 
X  +  e^x  +  e^ : 

e^  =  C3  -  e^d^  -  d^ 


^2  =  ^2  -  ^-^2  -  \ 
\-  c^-   e^dj_ 

To  obtain  a  set  of  formulas  for  b's,  d's  and  e's  from  which  the 
iterative  formulas  are  to  be  based,  the  idea  is  to  set  the  remainder 
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coefficients  r's  and  r's  to  0  and  eliminate  the  c's  using  the  above 
equations.  Thus  by  first  setting  r  and  r  to  0^  we  get 

d,  =  —  and  d^  =  . 

1   e^       2      e^ 

to  eliminate  c,  and  Cp,  use  the  first  group  of  formulas  above  to  obtadn 

5^  -   V     -   a^b|  -  hi  *   Jb^bj  -   a^bj^  -  e^d^ 
^2-  q ■ 

Similarly,   by  setting  r^    and  r_  to  0,    and  eliminating  the  c's, 

a^  -  b,  e,d^  -  b,  e^d. 


\           K          ^2  -  '^1^1^2  -  ^^2^1         , 
\    =  T^-  >     b„   =  -r— and 


1       d^e^>        2  d^e^ 

^2   =  ^6  -  ^2   -  ^2'      ^1  =   ^5  "  ^  ~  ^2^6  ^  ^2   -  ^  -  <^2^2- 
So  the  iterative  formulas  can  now  be   set  up  for  the  b's,    d's 
and  e ' s  in  terms  of  their  previous  values  and  of  the  original  coeffi- 
cients. 

For  the  case  of  linear  factors,    (x+x. ) ,   i  =  1,    . . .  n,   the 
fonnulas  are 

X     =  a     -  Ex. 
n         n  X 

X     T    =   [a     -,    -   a  Ex.    +  E(x. +x.x.)  ]/x 
n-1  n-1         n     i  '   i     i  J   -^     n 


X     o  =  [a     o  -   a     ^Ex. +  a   E(x,+x,x.)   -  E(x.+x.x,)  ]/x     ^x 
n-2       ■•   n-2         n-1     i      ni     ij  i     ij         n-1  n 
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n-3     n-3    n-2  i    n-1  '  i  x  y         n   i  i  j^ 


)    ^ 

+  i;(x.+  xrx.)]/x    ^  X   ^  X 

'  1  X  y^     n-2  n-1  n 


\  -   VV3  •••  \ 

where  the  notation  T.  xs   used  in  the  symmetric  function  sense,  it  being 

understood  that  in  the  formula  for  x  , 

P 


!C.  +  X.    X  . )  =    Z    (x.  +  X.    X. 

1     1     y    i,j=i   1    i     J 


r  =  a^  -  pb^  -  qb^  -  ^^ 


h.     Bairstow'  s  Method  for  Determining 
a  Quadratic  Factor  [6 J 

From  the  recursive  relations  among  the  coefficients  of  the 

n-1    .  ^ 
quotient   E  b^x^~  of  f (x)  by  the  quadratic  factor  x^  +  px  +  q, 

namely, 

\-k-i  =  Vk.i  -  pVk  -  ^Vk+i'  ^  =  1'  •••  "-2^  Vi  =  1'  \  =  ^> 

we  get  the  remainder  terms  with  coefficients 

\ 
s  =  a^  -  qb^       =  b_^  +  pbQ 

where  b^  and  b_^  are  obtained  by  taking  k  =  n-1,  n,  respectively. 

V/e  require  r  and  s,   i^ich  are   actually  functions  of  p  and  q, 
say  f-j_(p,q)    and  f2(p,q),    respectively,    to  go  to  0.      To  do  this, 
express  the  functions  f^  and  f^   as  Taylor  Series  in  two  variables, 
neglecting  all  terms  beyond  the  first  three,    and  equate  to  0: 
0  =  fi(p,q)    =  \{t^,\)   *  Ap  fip(Pj^,qj^)    +  Aq  f^qCpj^^qj^)      and 
0  =  f2(p,q)   =  i"2(Pk.\)   +  Ap  S^^{^^,c^^)    ^  Aq  f2q(Pj^,qi^)     where 
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(p,  ,q,  )  is  the  estimate  for  (pjq)  at  iteration  k,  Ap  =  p    -  p  , 

6f. 
^^  "  \+l  ~  \'   ^ij^Pk'\^  denotes  -^  evaluated  at  (Pj^jq^^)- 

Ap  and  Aq  may  be  solved  for  simiiltaneously  from  the  above  two 
equations  as  follows: 

Substituting  the  values  for  the  partials  of  f^   and  f^  with 
respect  to  p  and  q,  as  obtained  from  the  equations  for  r  and  s, 

6b  ,     6b„  6b  ,     6b 


0  =  b  ^  +  pb^  +  ^p(-r^)  +  p-r^   +  b.)  +  Aq(--li  +  p  -^) 
-1^0    ^  Sp  ^   '^Op    0'     ^  6q    ^  6q  ' 


Adding  (-p)  times  the  first  equation  to  the  second,  we  obtain 

%  *^1 

(1)  0  =  b^  .  Ap  ^  .  Aq  -^ 

6b  6b_ 

(2)  0  =  b_^  -  ApC-jZi  .  b^)  .  Aq  -^  . 

Now  to  obtain  the  partials  of  b  and  b_  ,  we  may  use  the  recursive 
relations 

6b  ,  ,  6b  ,      Sb  .  , 

n-k-1  n-k   „   n-k+1   ,       ^„  , 

— T. =  -  P  — 5 -  q  r -  b  ,    and 

6p  ^   Op    ^    oq      n-k 

^b  ,  .  6b  ,  6b  ,  , 

n-k-1  n-k  ^   n-k+1   ,       ^*,^„v, 

. =  -  p  — r -  q  5 -  b  ,  T  wnicn 

Oq  ^   Oq  ^    Oq      n-k+1 

are  obtained  from  the  recursive  formulas  of  the  b's  above.  However, 

for  computational  purposes,  it  is  worthy  of  note  that  if  the  quotient 

n-1    ,  ,  2 

y.     b  X    is  further  divided  by  x  +  px  +  q,  we  get  as  before  the 
1   ^ 
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recursive  relation  c  ,  ^=b  .^-pc  ,  _  -qc  ,  ^   for  the 
n-k-3    n-k-1  '^   n-k-2  ^  n-k-1 

coefficients  of  the  quotient.  Replacing  k  by  k-l^  and  comparing  with 

the  first  recursive  relation  for  the  partial  of  b's  (that  is,  with 

respect  to  p),  we  see  that 

6b 


n-k-1 


3^     "  Vk-2' 


similarly,  replacing  k  by  k-2,  and  comparing  with  the  second  (that  is, 
partials  with  respect  to  q) , 

6b 


n-k-1  ^ 
6q     ~  n-k-1, 


Thus   (1)    and   (2)    above  reduce  to 

(1)  b_^  +  ^P(-c_2   +  ^q)    +  Aq(-c_^)    =  0 

(2)  b^     +  Ap(-c_^)    +  Aq(-CQ)  =  0 

from  which  we  derive  the  formulas 


Ap  =  —, and 

"-1  "  "0^V^2^ 


Aq  =  — 


^-1^-1  ■*■  ^O^V^-2) 


c^l  ^  CQ(bQ  -  c_2) 


If  |Ap|  +|Aq|  <  e,  a  small  preassigned  number,  then  the  iter- 
as  converged  to  the  desired  p  and  q; 
and  q    =  Aq+q   and  repeat  the  procedure, 


ation  has  converged  to  the  desired  p  and  q;  otherwise,  let  p,  ,  =  Ap+q 

'     ^k+1    ^  ^k 
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To  recapitulate  the  procedure: 

(a)  Take   a  starting  value  for  p  as  well  as  for  q. 

(b)  Obtain  b     ^,    ...  b^.b  ^   using 

n-ti  0     -1 

^-k-1  =  \-k+i  -  P^n-k  -  ^\-k^r 

k  =  1,    . . .  n  vrtiere  b     =  0  and  b     ,    =  1. 
n  n-1 

(c)  Derive  c  i  ,  ...  c^,  c  t,  c  „  from  the  relation 

n-k  0   -1   -2 

^n-k-3  ^  Vk-1  -  P'=n-k-2  -  ^'^n-k-l' 
k  =  1,  ...  n-1. 

(d)  Use  C-,  c  ^,  c  ^  to  obtain  Ap  and  Aq  from 

Ap  =  -X and 

^-1  ^  °0^^0  -  ^-2^ 


,    ^-l"-l  "  ^O^V^-2^ 
Aq  =  — ^ 


"-1  "  "0^V"-2) 


(e)  Test  for  |Ap|  +  |Aq| 


If  <  e,  take  (p  ,q  )  to  be  the  desired  (p,q)  and  solve 

2 
the  equation  x  +  px  +  q  for  the  zeros.  Otherwise,  repeat  the 

iteration  using  p.  ,  =  Ap  +  p  and  q.  ^  =  Aq  +  q  ,  unless  the 

number  of  iterations  has  exceeded  the  upper  limit  set  at  the  start. 


5.  Graeffe's  Root  Squaring  Method  for  f(x) 
Having  Complex  Zeros  [12 J 

n+1         .  _, 

From  f(x)    =     E     a.x       ,    derive  g(x)    =  f  (v'x  )  -iriiose   zeros 
1       ^ 
are   squares  of  the  zeros  of  f (x) .     For  instance,   for  n  =  k, 

X     +  a,  x     +    .  .  .    +  a^    is  transformed  into  x     +   (-a,    +  2a^)x     + 
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2     2  2  2 

(2a.^    -  2a.  a     +  a  )x     +   (2a^a     -  a^)x  +  a^  so  that  if 

"^^         1-1 
g(x)   =     E     b.x       ,   we  have  the  formulas 

1       ^ 

\  =  -%  ^  2a3 

b3  =  a^  -  2a^a2  +  2a^ 
^2   =   -4  "   2a3ai 

^   =   4. 

If  the  zeros  of  f(x)   differ  in  absolute  value,   the   zeros  of 
g(x)  must  differ  by  a  greater  amount.     Thus  by  doing  this  transforma- 
tion m  times,    an  equation  may  be  foiind  vdth  roots  equal  to  the  2mth 
power  of  the  original  and  which  roots  differ  in  absolute  value  as 
greatly  as  may  be  desired.     In  this  case,   the  coefficient  of  x 
will  nearly  be  the   square  of  the  absolute  value  of  the  zeros  with 
largest  magnitude,    since  if  a  +  bi  are  those   zeros,   the  S^r.r. 

2  2 

is  dominated  by  a     +  b    .     The  greater  m  is,    the  more  nearly  true  this 

is.     The  root-squaring  process  is  terminated  when  the  coefficient  of 

X         in  the  last  derived  equation  is  nearly  equal  that  in  the  previous 

step. 

Having  found  in  this  manner  the  Umth  power  of  the  largest 

absolute  value  among  the  zeros  of  the  original,  the  square  of  such 

absolute  value  may  be  obtained  by  logarithm.  The  zeros  can  then 

be  determined  from  this  by  using  relation  among  the  zeros. 
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To  illustrate  for  n  =  U,  if  the  zeros  are  c  +  d.i,  J  =  1,2, 

P      P     P      P  0  0  0 

then  a^  =  (c^  +  ^{^^^  ""  '^2^   ^°  ^'^^''^  since  |r^|   =  c^  +  cL  is  known, 

the  square  of  the  absolute  values  of  the  second  pair  of  zeros, 

I   |2    2    2 

'^2   "^  ^2  "*"  *^2  *^^"  ^®  determined.  Then  using  the  relation 

\   =  -2(c^  +  c^)  and  a^  =  -2(|r^|  c^  +  Ir^l  c^),  c^  and  c^  can 
be  solved  for  simultaneously,  and  consequently  from  the  values  of 

|r  I   and  |r_|  ,  d^  and  d  can  be  determined. 


6.     Hitchcock's  GCD  Method  for  f(x) 
Having  All  Zeros  Complex  [14,5] 

First  assume  the  coefficients  a.    of  f(x)   real. 

"■"^         i-1  2 

By  dividing  f(x)    =     E     a  x         by  x     +  q,   we  obtain  the  identity 

1   ^ 

f(x)  =  [x""^  +  a  x"-^  +  (a  ,  -  q)x"-^  +  ...  ][x^+q]  + 
n        n-1 

r  2  m-l,    r  2      —  m, 

^^2  ~   ^\^  ^   ^b°^     -  .  •  •  1  a^q   ]x  +  [a^  -  a^q  +  a^q  -  .  . .  +  q  ] 

where  m  =  n/2  and  the  +  in  the  last  term  of  each  of  the  last  two  coeffi- 
cient expressions  is  determined  such  that  the  signs  are  alternating. 

If  +  J^   were  a  pair  of  zeros  of  f (x) ,  then  the  last  two  coef- 
ficients expressions  would  be  identically  0.  Hence  we  derive  the  two 
test  equations  for  a  pair  of  pure  imaginary  zeros,  namely, 

/-iN  2  m-1   _ 

(1)  a^  -  a^^q  +  a^q  -  ...  +  a^q    =  0 

I  r^\  2      —  m       - 

(2)  a^  -  a^q  +  a^q  -  ...  +  q      =0 
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So  given  f(x),  if  the  above  polynomials  (1)  and  (2)  have  a  common 
factor,  say  q-q.,  ,  then  we  shall  have  determined  a  pair  of  zeros  of 
f(x),  namely,  i  v'qT  i.   The  GCD  algorithm  can  be  employed  for  this 
purpose.  However,  in  most  cases  the  common  factor  is  1.  So  the  idea 
now  is  to  use  Horner's  root  diminishing  process  and  move  a  pair  of 
zeros  of  f(x)  until  they  are  transformed  to  pure  imaginary  ones. 

The  process  then  boils  do;m  to  using  the  GCD  algorithm  on  the 
test  equations.  If  the  last  real  number  remainder  in  the  algorithm 
is  positive  (or  negative),  try  translation  of  variable  by  a  pegitive 
(or  positive)  number  h,  and  apply  the  GCD  algorithm  again  on  the  new 
test  equations.  If  the  remainder  changes  sign,  then  we  know  there  is 
a  pair  of  zeros  whose  real  part  lies  between  0  and  +h.  So  try  a  num- 
ber between  0  and  h  for  a  translating  factor,  and  repeat  the  process. 
Continuing  this  way,  we  can  thus  narrow  down  the  interval  containing 
the  desired  real  part  by  translations  till  the  remainder  in  the  GCD 
algorithm  approximates  0.  This  way,  we  also  find  at  the  same  time  the 
imaginary  part  of  that  particular  pair  of  zeros  from  the  common  factor 
of  the  two  test  polynomials.  We  may  also  employ  interpolation  to 
estimate  the  amount  of  translation  necessary,  after  we  have  tried  a 
couple  of  values  for  h.  Thus  if  translation  of  variable  by  h  gives 
a  common  factor  x-q  to  the  test  polynomials,  the  desired  zeros  are 
then  +h  +  Vq  i. 

For  a . ' s  complex,  the  procedure  is  practically  the  same  except 

that  there  is  a  different  set  of  test  equations  for  a  pure  imaginary 

zero  qi.  If  f(x)  =  x"  +  (a  +  bi)x""^  +  ...  (a  +  hi),  then  the 

n    n  -1-1 
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remainder  in  dividing  f(x)  by  x-qi  is  f(qi),  which  is  what  we  need  to 
test  whether  0  or  not.   By  separating  real  and  imaginary  parts  of  fCqi) 
the  two  test  equations  then  are 

(1)  a^  -  b^q  -  a^q  +  b^q  +  a^q  -  b^q^  -  a^q  . , .  =  0 

(2)  b^  +  a^q  -  b^q  -  a^q  +  b^q  +  a^q^  -  b^q  -...  =0 

where  it  will  be  noted  that  the  signs  are  2  plus  and  2  minus  alter- 
natively except  at  the  start  of  (1) . 
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